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'Electromagnetic mass" where gravitational mass and other physical quantities orig- 



Q I inate from the electromagnetic field alone has a century long distinguished history. 

In the introductory chapter we have divided this history into three broad categories 

" I 

^ ' - classical, quantum mechanical and general relativistic. Each of the categories has 

^ been described at a length to get the detailed picture of the physical background. 

Recent developments on Repulsive Electromagnetic Mass Models are of special in- 
terest in this introductory part of the thesis. In this context we have also stated 
motivation of our work. In the subsequent chapters we have presented our results 
and their physical significances. It is concluded that the electromagnetic mass mod- 
els which are the sources of purely electromagnetic origin "have not only heuristic 
flavor associated with the conjecture of Lorentz but even a physics having uncon- 
ventional yet novel features characterizing their own contributions independent of 
the rest of the physics" . 
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SYNOPSIS 



The investigations carried out in the thesis "ELECTROMAGNETIC MASS 
MODELS IN GENERAL THEORY OF RELATIVITY" form seven chapters in- 
cluding the introductory and concluding ones. 

Introduction 

The study of "electromagnetic mass" where gravitational mass and other physical 
quantities originate from the electromagnetic field alone has a century long distin- 
guished history. In the introductory chapter we have divided this history into three 
broad categories - classical, quantum mechanical and general relativistic. Each of 
the categories has been described at a length to get the detailed picture of the physi- 
cal background. In the classical part starting from the Lorentz's Theory of Electrons 
it includes Thomson's Concept of Electromagnetic Mass and Abraham's Model for 
the Electrons. We have also described the drawbacks of Lorentz's Model for the Elec- 
trons. In this connection Poincare's Theory of Electrons has been put forward to 
eliminate the discrepancy in Lorentz's Model. Einstein's Special Relativistic Model 
of Electrons and various other Models including the General Relativistic one and 
Quantum Electron Models have been included in other parts of the history. Recent 
developments on Repulsive Electromagnetic Mass Models are of special interest in 
this introductory part of the thesis. In this context we have also stated motivation 
of our work. 

CHAPTER II 

Relativistic Electromagnetic Mass Models with Cosmological Variable A 

in Spherically Symmetric Anisotropic Source 

In the chapter II a class of exact solutions for the Einstein-Maxwell field equa- 
tions, which are obtained by assuming the erstwhile cosmological constant A to be 
a space-variable scalar, viz., A = A(r). The source considered here is static, spheri- 
cally symmetric and anisotropic charged fluid. The solutions obtained are matched 
continuously to the exterior Reissner- Nordstrom solution and each of the four solu- 
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tions represents an electromagnetic mass model. 
CHAPTER III 

Classical Electron Model with Negative Energy Density in Einstein-Ccirtan 
Theory of Gravitation 

The experimental results regarding the maximum limit of the radius of the electron 
~ 10~^^ cm and a few of the theoretical works readily suggest that there might be 
some negative energy density regions within the particle as per General Theory of 
Relativity. It is argued in the chapter III of the present investigation that such a 
negative energy density can also be obtained with a better physical interpretation 
in the framework of the Einstein- Cartan theory. 

CHAPTER IV 

Energy Density in General Relativity: a Possible Role for Cosmological 
Constant 

We consider a static spherically symmetric charged anisotropic fluid source of ra- 
dius ~ 10~^^ cm by introducing a variable A dependent on the radial coordinate r 
under general relativity in the chapter IV. Prom the solution sets a possible role of 
the cosmological constant is investigated which indicates the dependence of energy 
density on it. 

CHAPTER V 

Relativistic Electromagnetic Mass Models: Charged Dust Distribution 
in Higher Dimensions 

Electromagnetic mass models are proved to exist in higher dimensional theory of 
general relativity corresponding to charged dust distribution. In the chapter V, 
along with the general proof, a specific example is also cited as a supporting candi- 
date. 
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CHAPTER VI 

Relativistic Anisotropic Chcirged Fluid Spheres with Veirying Cosmolog- 
ical Constant 

Static spherically symmetric anisotropic source has been studied for the Einstein- 
Maxwell field equations assuming the erstwhile cosmological constant A to be a 
space- variable scalar, viz., A = A(r). The solutions thus obtained are shown to 
be electromagnetic in origin in the sense that all the physical parameters including 
the gravitational mass originate from the electromagnetic field alone. Moreover, to 
construct the models it is also shown that the generally used pure charge condition, 
viz., p + Pr = is not always required for constructing electromagnetic mass models. 
This is the main theme of the chapter VI. 

CHAPTER VII 
Conclusions 

Electromagnetic mass models which are the sources of purely electromagnetic origin 
"have not only heuristic fiavor associated with the conjecture of Lorentz but even a 
physics having unconventional yet novel features characterizing their own contribu- 
tions independent of the rest of the physics" (Tiwari 2001). This is, as Tiwari (2001) 
guess "may be due to the subtle nature of the mass of the source (being dependent 
on the electromagnetic field alone)" . Therefore, in our whole attempt we have tried 
to explore "the subtle nature of the mass of the source" . However, to do this under 
the general relativistic framework, we have considered Einstein field equations in its 
general form, i.e., with cosmological constant A which also acts as a source term to 
the energy-momentum tensor. If we consider that A has a variable structure which 
is dependent on the radial coordinate of the spherical distribution, viz., A = A(r) 
then it can be shown that A is related to pressure and matter energy density. Hence 
it contributes to the effective gravitational mass of the system. 
It is seen that equation of state has an important role in connection to electromag- 
netic mass model. Therefore, at first we have obtained electromagnetic mass model 
under the condition p -\- p — 0. However, later on it is shown that electromagnetic 
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mass model can also be obtained by using more general condition p + p 7^ 0. 
The model considered in our work, in general, corresponds to a charged sphere with 
cosmological parameter in such a way that it does not vanish at the boundary. The 
idea behind is that the cosmological parameter is related to the zero point vacuum 
energy it should have some finite non-zero value even at the surface of the bounding 
system. For this type of spherical system we can have a class of solutions related to 
charged as well as neutral configurations. 

It can be shown that these models have positive energy densities everywhere. Their 
corresponding radii are always much larger than 10~^^ cm. Furthermore, as the 
radii of these models shrink to zero, their total gravitational mass becomes infinite. 
It have been shown by Bonnor and Coopcrstock (1989) that an electron must have 
a negative energy distribution (at least for some values of the radial coordinate). 
In this connection we have shown that the cosmological parameter A has a definite 
role on the energy density of micro particle, like electron. At an early epoch of 
the universe when the numerical value of negative A was higher than that of the 
energy density p, the later quantity became a positive one. In the case of decreasing 
negative value of A there was a smooth crossover from positive energy density to a 
negative energy density. 

It is suggested by Bonnor and Cooperstock (1989) and Herrera and Varela (1994) 
that spin and magnetic moment can be introduced to electron through the Kerr- 
Newman metric. But, it has been seen that the Kerr-Newman metric cannot be 
valid for distance scales of the radius of a subatomic particle. We, therefore, tackled 
the problem in the frame work of Einstein-Cartan theory where torsion and spin are 
inherently present. In this case, the only way is to take the spin to be the 'intrinsic 
angular momentum' that is the spin of quantum mechanical origin. In our work 
considering the spins of all the individual particles are assumed to be oriented along 
the radial axis of the spherical systems we have obtained some interesting solutions 
with physical validity. 

Another important point we would like to mention here that in all the previous inves- 
tigations we have studied electromagnetic mass models in 4-dimensional Einstein- 
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Maxwell space-times only. Therefore, one can ask whether electromagnetic mass 
models also can exist in higher dimensional theory of General Relativity. We have 
presented a model which corresponds to spherically symmetric gravitational sources 
of purely electromagnetic origin in the space-time of (n -|- 2) dimensional theory of 
general relativity. 

We have also taken up the problem of anisotropic fluid sphere as studied earlier in 
a different view point. By expressing A in terms of electric field strength E we have 
explored some possibilities to construct electromagnetic mass models using the con- 
straint p + p ^ We would like to mention here that unhke the solutions of Gr0n 

(1986a,b) and Ponce de Leon (1987a,b) in the present investigation, in general, the 
electric field (and hence the cosmological constant) does not vanish at the bound- 
ary. However, it is shown that the class of solutions obtained here are related to 
charged as well as neutral systems of Gr0n (1986a,b) and Ponce de Leon (1987a,b) 
depending on the values of the parameter N . 
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Chapter 1 
Introduction 



"...the world stands before us as a great, eternal riddle. 
- Albert Einstein 



The first elementary particle was called corpuscle by J. J. Thomson (1881), later 
named as electron - the word proposed by G. Johnstone Stoney in 1891 as a unit of 
electronic charge. There were two schools of thought, one school favoring atomistic 
world view, and the other believing in the continuum. The debate continued for 
several years in absence of sufficient experimental findings. It is in this respect that 
J. J. Thomson's work turns out to be epoch making, and the credit for the discovery 
of electron belongs to him. He succeeded in a determination of the charge to mass 
ratio and the elementary charge e. The physical world and all its phenomena in 
principle can be reduced to the problem of the interaction of the elementary parti- 
cles. The four types of fundamental interacting forces are labeled as electromagnetic, 
gravitational, strong and weak. Of these four interactions the electromagnetic one 
thoroughly investigated and has been best understood in connection with the elec- 
trically charged particles, especially, the electrons. 

The study of electromagnetic mass with a century long distinguished history can be 
divided into three broad categories - classical, quantum mechanical and general rel- 
ativistic. Here we shall first follow the tradition established by the classical electron 
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theory. Actually, there are two reasons why the classical electron is studied: 

a) There is no quantum mechanical model for the electron; 

b) It is quite natural to complete the classical electromagnetic mass theory first 
before making the transition to a quantum description of the inertial properties of 
an electron. 

In the classical framework Lorentz tried to tackle the problem of the electrodynam- 
ics of moving bodies. To get an overview of the problem starting from Lorentz, we 
shall first provide a brief historical account of the work done by different authors in 
a sequence. 

1.1 A Brief Historical Background of the Theories 
About the Structure of the Electrons 

1.1.1 Lorentz 's Theory of Electrons 

Lorentz's (1892) apparent motivation was to solve the null result of the Michelson- 
Morley experiment keeping the existence of ether as it is. He wanted to represent 
an electromagnetic world view in comparison to the Maxwellian electromagnetic 
theory. Based on this philosophy he developed his Theory of Electron in 1892. 
The basic assumption made by him is that microscopic charged particles or ions 
in motion through absolutely resting ether were the source of the electromagnetic 
disturbances. He further assumed that ether permeates the electrical particles and 
that electrical particles were perfectly rigid bodies. He considered an electric field 

— * — * 

vector E and a magnetic field vector H in this absolute ether frame and obtained 



where a — > electric charge density, u — > velocity of electric particles and c — > velocity 
of light quanta, photons with velocity c = 2.99 x lO^^cm/sec in vacua. 




(1.1) 
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He also considered the Maxwell equations 

V.^ = 47rp, (1.2) 
V.B = 0, (1.3) 
Vx^=-if (1.4) 

and 

where E, B, J and c respectively are the electric field, magnetic field, current density 
and velocity of light in vacua. 

Now, assuming that the Lorentz- Maxwell equations (1.1-1.5) were valid in the rest 
frame of the ether, Lorentz considered the Galilean transformations 

X — X — vt, y — y, z — z, t — t (1-6) 

for a reference frame S moving with an uniform velocity v in the direction of x 
with respect to the frame S. To retain the form of the wave equation describing an 
electromagnetic disturbance in both the reference frames S and S Lorentz introduced 
new variables x', y', z' and t' connected with the equations 

x' = X, y' = y, z' = z, t' = t — -x (1.7) 

(with upto first order of v/c). 

Applying this theory Lorentz was able to explain the electrodynamics of moving bod- 
ies and various optical phenomena like the propagation of light in dielectrics at rest 
as well as in moving media and hence the Fresnel's dragging effect. The first fruit- 
ful idea for explaining the null result of Michelson-Morley came from Lorentz. He 
suggested that material bodies contract when they are moving, and the shortening 
is only in the direction of the motion. He proposed that the length of the interfer- 
ometer arm parallel to the direction of motion is shortened by a factor (1 — v^/2c^). 
This is known as Lorentz-Fitz Gerald contraction hypothesis. The more simplified 
and generalized theory of Lorentz (1895) on electrodynamics which involved the 
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concept of local time in contrast to the universal true time was put forward by him 
to ensure that the form of Maxwell equations for charge-free-space remains the same 
in a moving frame (at least upto first order inv/c). Lorentz assumed that there are 
several electrons in each atom which are elastically bound to an equilibrium position 
and thus are able to perform harmonic vibrations with given frequencies. In electric 
conductors additional electrons were assumed to move freely. With these funda- 
mental theoretical tools it was possible to explain a great number of phenomena - 
the absorption, scattering and refraction of light by matter, the Zeeman effect and 
many more. During 1900 to 1903 Lorentz conjectured that a part of the electron 
mass might be of electromagnetic origin. However, for the sake of history we should 
mention that even before Lorentz there were other notable scientists, who had ex- 
pressed the idea of electromagnetic mass in their works. 

1.1.2 Thomson's Concept of Electromagnetic Mass 

While studying the interaction of charged particles Thomson (1881) found that the 
kinetic energy of a charged sphere increases by its motion through a medium of finite 
specific inductive capacity. He pointed out that the increase in the kinetic energy 
was due to the self induced magnetic field of the charged sphere and calculated the 
total mass to be 

m = mo + II (1.8) 

with 11 = ae^ jac? where mo is the mass of the charged sphere, \x is the increased 
mass, e the electric charge, a the radius of the sphere and a is a numerical factor 
of order unity. Thus he came to the conclusion that "the effect of electrification 
is the same as if the mass of the sphere were increased This increased mass 
which clearly is of electromagnetic origin was termed as electromagnetic mass. The 
works of Thomson on electromagnetic mass were improved upon by O. Heaviside 
(1889) showing that the mass of a uniformly moving charged body varied with ve- 
locity. G. F. C. Searle (1897) extended the work of Heaviside as that the energy of a 
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charged body and therefore its mass increases with velocity. Later on W. Kaufmann 
(1901 a,b,c; 1902 a,b) through a series of experiments on beta rays estabhshed the 
dependence of the electron mass on velocity. He (1901b) showed that one-third of 
the fast moving electron mass was of electromagnetic origin. Hence it was proba- 
bly an easy task for M. Abraham (1902) to speculate that the electron mass was 
completely of electromagnetic origin. By reanalyzing his experimental data Kauf- 
mann (1902b) found that the mass of the electron is purely of electromagnetic origin. 

1.1.3 Abraham's Model for the Electrons 

M. Abraham (1902, 1905) proposed the first field theoretical model for the electron 
on the basis of Lorentz field equation (1.1) and the Maxwell field equations (1.2 
-1.5). He assumed the electron to be a rigid sphere with uniform surface charge 
distribution. Further, he considered the electron charge density as a fundamental 

quantity. From Kaufmann's experimental work he ah^cady knew about the velocity 
dependence of electron mass. So, he undertook a detailed study of the dynamics of 
electron including a scheme to derive the electron mass entirely from its self-field. 
Proposing the equation of motion for electron as an analogy for Newtonian equations 
of motion and considering that the total force acting on the electron should always 
vanish, Abraham ultimately got two types of electron masses, one is the longitudinal 
mass (my) and the other is transverse mass {m±) as 




(1.9) 



and 



(l + /3')/n 



1+/9 



l-(3 



) 



-2(3 



(1.10) 



For low velocities and in the limit /? = 0, the electron mass was 




(1.11) 
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He also established the electron radius to be of the order of 10~^^ cm. Finally, he 
claimed that the transverse mass (1-10) was in agreement with Kaufmann's experi- 
mental data to a good approximation. 



1.1.4 Lorentz's Model for the Electrons 

In 1903 Lorentz extended his electromagnetic theory of 1895. Here he included the 
equations of motion of free electron together with a review of Abraham's model of 
the electrons. In the further development of the theory Lorentz speculated that the 
spherical electron would experience an ellipsoidal change in its shape while it is in 
motion and obtained the transverse mass (the relativistic mass) for electron as 

^^=67rac^(l'-/3^)V2- ^^'^^^ 
It is then a straight forward way to obtain the relativistic electromagnetic mass of 
the field of the electron as 

which is not the same as non-relativistic electromagnetic mass 

™.,.c = 5 (^) . (1.14) 

Combining these two one can write 

Uelec = ^melecC^. (1.15) 

Certainly this relation unifies gravitation with electromagnctism and thus through 
Lorentz the concept of electromagnetic mass was born. Lorentz then proposed a 
model of an electron as an extended body consisting of only pure charge and no 
matter and the charge is uniformly distributed on a spherical shell. In his work 
Lorentz dealt only with the inertial aspect of mass. Lorentz was certainly aware 
of the Newtonian gravitational aspect of mass, but he probably disregarded grav- 
itational effects because Newtonian gravitational forces are smaller than electrical 
forces by many orders of magnitude. 
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1.1.5 Drawbacks of Lorentz's Model 



A discrepancy lies between the two formulas (equations 1.13 and 1.14) for the elec- 
tromagnetic mass. This problem implies that the relationship between momentum 
and velocity for a particle in Newtonian mechanics differs from that for the com- 
pletely electromagnetic electron. This defect can be corrected by merging this theory 
with special relativity. The anomaly of factor 4/3 disappears since it is incompatible 
with the rclativistic transformation properties. For the finite electron this was first 
pointed out by Fermi (1922), but his work did not receive its recognition till 1965. 
For point electrons the removal of this factor was later rediscovered. The inertia 
and mass of the classical electron originate from the unbalanced mutual repulsion 
of the volume elements of the charge caused by the distorted electric field of an 
accelerating electron. However, it is not clear what keeps the electron stable since 
the Lorentz's model of the electron describes its charge as uniformly distributed on 
a spherical shell, which means that its volume elements tend to blow up by repelling 
one another. This difficulty can be removed by eliminating the electron structure 
and assuming the particle to be a point particle. This indeed produces a new dif- 
ficulty, i.e., when the radius shrinks to zero, the electron's mass becomes infinitely 
large. This is the famous self energy problem. It exists in the classical theory as well 
as in the quantum theory of the electron. Its satisfactory solution is not yet known. 
Thus the electron is then considered to have a finite extension. The difficulty with 
the Lorentz's model was that it had no mechanism to overcome the electrostatic 
repulsion of the charge, so that the body was unstable. 

1.1.6 Poincare's Theory of Electrons 

Poincare (1905, 1906), with the aim of overcoming the instability and inconsis- 
tency of Abraham's model with respect to the special relativistic Lorentz transfor- 
mations proposed that an attractive force of cohesive type and consequently non- 
electromagnetic in nature can be added so as to just balance the stresses and estab- 
lish stability. He first considered the coordinate transformations in the proper form 
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and called these as Lorentz transformations which are 



x' = ^{x — vt), y' = y, z' = z, t' = ^{t — vx/(?). 



(1.16) 



Then he proposed a new Lagrangian which was composed of two parts 



L — L field + L stress 



(1.17) 



where Lfieid the Lagrangian for the electron's self electromagnetic fields, L stress ~^ 
the Lagrangian related with the so-called Poincare stress and Mq"* electromag- 
netic energy of the spherical body. 

This new Lagrangian not only solved the instability problem but also removed the 
discrepancy in the calculation of mass of the electron. In fact, one can show that in 
a relativistic theory in which the electron self energy is not infinite, the self stress 
will vanish and the particle will be stable. 

1.1.7 Einstein's Special Relativistic Model of Electrons 

Einstein's special theory of relativity is based on two unique postulates, that is, the 
principle of relativity and the constancy of velocity of light in vacuum. Using these 
postulates he could derive the Lorentz transformations and explain the length con- 
traction and time dilation as a kinematical consequence of these transformations. 

He further showed the covariance of Maxwell-Lorentz electromagnetic field equa- 
tions under Lorentz transformation and that the Lorentz force equation (1.1) is a 
consequence of the principle of relativity. The longitudinal and transverse masses 
of the electron as given by Einstein (1905) are 



mil = 7^mo 




and 



m_L = 7mo 



(1.19) 
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where tuq is the rest mass of the electron. These results are exactly identical to 
those of Lorentz. Of course, the dependence of mass on velocity, as given by equa- 
tion (1.19) does not mean that only the electromagnetic mass of an electron is a 
privileged mass to vary with velocity; rather it is just like any other kind of mass. 

1.1.8 Various Other Type of Models Including General Rel- 
ativistic Electron Models 

In 1912 G. Mie (1912 a,b) tried to build charged particle models based on electro- 
magnetic fields alone, so that the mass of the charged particle hke electron could be 
completely of electromagnetic origin, and suggested a modification of the Maxwell- 
Lorentz field equations. He assumed that the complete electromagnetic field is 
determined by ten universal quantities which are functions of the four-potentials 
and the Maxwell tensors F*-'. But this unitary field theory ultimately failed. 
The other workers on the theory of unitary field, with different view points, were 
M. Born and L. Infeld (1934), B. Hoffmann (1935 a,b) and F. Bopp (1940, 1943). 
The search for a solution to the problem of the electron structure were made by H. 
Weyl (1918a,b; 1919), T. Kaluza (1921) and O. Klein (1926a,b; 1928) in the realm 
of unified field theory where gravitational and electromagnetic fields have been uni- 
fied into a single theory. In 1919, the first general relativistic approach towards an 
electromagnetic mass was put forward by Einstein. To overcome the drawbacks of 
Mie's theory Einstein proposed a model where gravitational forces would provide 
the necessary stability to the electron and the contribution to the mass would also 
come from it. The discovery of quantum mechanics in 1925 - 1926 called for an ex- 
tension of the classical theory of the electron to the atomic and subatomic domain. 
A very important and unexpected property of electrons was found almost simul- 
taneously with the establishment of this new mechanics. It was discovered that 
the electron has an intrinsic angular momentum, a spin and a magnetic moment 
associated with it. Essential progress towards a relativistic quantum mechanical 
description of electrons was made by Dirac (1928). It was Dirac who was able to 
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devise a wave equation for the electron which fulfilled the relativistic requirements. 
Actually quantum mechanics treats electron as a point-like charged particle with 
spin and hence extended electron could not be accommodated within it. Thus it 
seems that for a better description of the electron structure with a spin in the Gen- 
eral Relativistic framwork it is already seen that the Einstein-Cartan-Maxwell or 
Einstein-Maxwell-Dirac space-times is preferable rather than the Einstein-Maxwell 
space-times (Ray and Bhadra 2004a). In the various classical point charge theories 
(Mehra 1973) the electron is treated as a point charge having a pure charge without 
any structure. But these theories cannot overcome the self energy problem of the 
electron which becomes infinite at its location. This infinite self energy problem can 
be solved by considering an extended charge distribution for the electron. But these 
theories have no satisfactory quantum versions. In the framework of Einstein's gen- 
eral theory of relativity (which was proposed in 1916) a lot of work has been carried 
out by different authors on charged body. Some models which were developed to 
study the structure of electron are due to Kyle and Martin (1967), Cohen and Cohen 
(1969) and Bayhn and Eimerl (1972). Katz and Horwitz (1971) and Lopez (1984) 
have developed classical extended electron models from the general relativistic point 
of view. In all these models an electron has been considered to be a microscopic 
sphere of charged perfect fluid or as a spherical shell of matter embodied by charge. 



1.2 A Short Account of the Recent Developments 
in the Electromagnetic Mass Models 

So long we have seen in the theories that in the Special Relativity and even in the 

General Relativity the mass of the electron lias been considered to contain two parts, 
that is, the non-electromagnetic part and the electromagnetic part. The exceptions 
are Lorentz's and Abraham's theories which have independently considered the mass 
of electron to be completely of electromagnetic origin. In the following section 
we shall give a brief description of the above idea showing how the conjectures of 
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Lorentz- Abraham may be revived on the ground of the General Relativity. 

1.2.1 Repulsive Electromagnetic Mass Models: Electron Type 

Repulsive gravitation is produced by the negative mass of the polarized vacuum. The 
vacuum fluid obeying an equation of state p = —p was taken by most of the workers 
for the construction of electromagnetic mass model. By considering the relation 
between the metric coefficients i.e., gooQu — ~1 (which for both the Schwarzschild 
and the Reissner- Nordstrom matrices equivalent to the relation 1/ + X — 0) to he 
valid inside a charged perfect fluid distribution, it is shown by Tiwari et al. (1984) 
that the mass energy density and the pressure of the distribution are of electromag- 
netic origin. In the absence of charge, however, there exists no interior solutions. A 
particular solution which conflrms the same and matches smoothly with the exte- 
rior Reissner-Nordstrom was obtained by them. This solution represents a charged 
particle whose mass is entirely of electromagnetic origin. The pressure being neg- 
ative here the model is under tension and hence the source is of repulsive nature. 
In the approach taken by Gautreau (1985), following Tiwari et al. (1984), the elec- 
tron's mass is associated with the Schwarzschild gravitational mass given by general 
relativity and not with the incrtial mass used by Lorentz (1904). In this case the 
Schwarzschild mass of an extended charged body as seen at infinity arises from the 
charge as well as the matter possessed by it. Here the field equations for a Lorentz 
type pure charge extended electron are obtained by setting the matter terms equal 
to zero in the field equations for a spherically symmetric charged perfect fluid. An 
explicit solution to the pure charge fleld equations are examined by Gautreau (1985). 
Lopez (1984) proposed a classical model of the spinning electron in which the particle 
is the source of the Kerr-Newman fleld. The electron here is regarded as a charged 
rotating shell with surface tension. The phenomenon of repulsive vacuum gravita- 
tion proved to be of importance in cosmology with appearance of the inflationary 
universe models. Gr0n (1985) pointed out the possibility that repulsive gravitation 
may be of importance also in connection with elementary particle models. This pos- 
sibility was realized by Tiwari et al. (1984) and by Lopez (1984). Poincare stresses 
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were explained by them as being due to vacuum polarization in connection with a 
recently presented class of electromagnetic mass models in general relativity. The 
gravitational blue shift of light is explained as being due to repulsive gravitation 
produced by the negative gravitational mass of the polarized vacuum. Gr0n (1985) 
pointed out that the electron model of Lopez (1984), which includes spin, and which 
is a source of the Kerr-Newman field gives rise to repulsive gravitation. 
Assuming an implicit relation among the unknown physical parameters viz., the 
pressure 'p, the charge density a and the electromagnetic potential 0, it has been 
shown by Tiwari et al. (1986) that <^ satisfies the well known Lane-Emden equa- 
tion. Electromagnetic mass models corresponding to the exact solutions of the 
Lane-Emden equation were obtained by them. The radii of some of the models were 
compared by them with the "classical electron radius". Lopez (1986) analyzed the 
stability of a classical ellipsoidal electron model. The model was found to be stable 
under oscillations which change the size of the ellipsoid without altering its shape. 
It is further shown by Lopez (1986) that angular momentum conservation does not 
allow the existence of other oscillation models. Ponce de Leon (1987a) investigated 
the relation g'oofl'ii = — 1 in the case in which the interior is filled with imperfect fluid. 
He found that the core of such a distribution is gravitationally repulsive provided 
the energy density is positive. Ponce de Leon (1988) also investigated the different 
aspect of the phenomenon of gravitational repulsion in static sources of the Reissner- 
Nordstrom field. He found that in the case of perfect fluid spheres there exists a 
close relation between the gravitational repulsion and the Weyl curvature tensor. 
Ponce de Leon (1988) proved that the static source of the Reissner-Nordstrom field 
gives rise to gravitational repulsion only if the pure gravitational field energy inside 
the sphere is negative. It is also proved that although the gravitational repulsion 
always takes place in the interior of a charged perfect fluid sphere when its radius 
is less than the classical electron radius, this is not necessarily so either in the case 
of anisotropic charged spheres or if the net charge of the body is concentrated at its 
boundary only. He further found that the charge contributes negatively to the effec- 
tive gravitational mass, in the sense that an increase in the charge causes a decrease 
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in gravitational mass. He explained the gravitational repulsion as being due to this 
negative contribution rather than the strain of vacuum because of vacuum polariza- 
tion. Bonnor and Cooperstock (1989) found, by modelling the electron as a charged 
sphere obeying Einstein-Maxwell theory, that it must contain some negative rest 
mass. The total gravitational mass within this sphere is negative which is one of the 
assumptions made in singularity theorems of general relativity. Lopez (1992) con- 
structed a classical model of the spinning electron in general relativity consisting of 
a rotating charge distribution with Poincare stresses. Obviously he obtained a class 
of interior solutions of the Kerr-Newman field. The negative pressures or tensions 
obtained here are identified with the cohesive forces introduced by Poincare (1905, 
1906) to stabilize the Lorcntz electron model. They arc shown by Lopez (1992) to 
be the source of a negative gravitational mass density and thereby of the violation 
of the energy conditions inside the electrons. Herrera and Varela (1994) pointed 
out the role played by the negative rest mass as mentioned in the work of Bonnor 
and Cooperstock (1989). Here the electron is modeled as a spherically symmetric 
charged distribution of matter deprived of spin and magnetic moment. Since the 
electrostatic energy of a point charge is infinite, the only way to produce a finite 
total mass is the presence of an infinite amount of negative energy at the center of 
symmetry. They (Herrera and Varela 1994), by analyzing some extended electron 
models, showed that negative energy distributions result from the requirement that 
the total mass of these models remains constant in the limit of a point particle. Ti- 
wari and Ray (1996) dealt with a model which is the charged generalization of static 
dust sphere in Einstein- Cartan theory. They obtained a set of solutions with torsion 
and spin which represents an electromagnetic mass model. Blinder (2001) proposed 
a model for the classical electron as a point charge with finite electromagnetic self 
energy. Modified form of the Reissner- Nordstrom and Kerr-Newman solutions of the 
electromagnetic equations were derived. Moreover, the self interaction of a charged 
particles with its own electromagnetic field was shown to be equivalent to its reaction 
to the vacuum polarization. 
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1.2.2 Repulsive Electromagnetic Mass Models: Stellar Type 

The work of Ray and Das (2002) which is concerned with the charged analogue 
of Bayin's work (1978) related to Tolman's type and presents astrophysically in- 
teresting aspects of stellar structure. However, in a static spherically symmetric 
Einstein-Maxwell space-time this class of astrophysical solution found out by Pant 
and Sah (1979) and Ray and Das (2002) has been revisited in connection with the 
phenomenological relationship between the gravitational and electromagnetic fields 
(Ray and Das 2004). Considering Riccati equation with known value of charge q for 
the total charge on the sphere in the following form 

q{a) = KoJ" (1.20) 

they have shown in one of the cases that the gravitational mass for n = 1 can be 

given as 



It is thus qualitatively shown that the charged relativistic stars of Tolman (1939) 
and Bayin (1978) type are of purely electromagnetic origin. Obviously, the exis- 
tence of this type of astrophysical solutions is a probable support to the extension 
of Lorentz's conjecture that electron-hke extended charged particle possesses only 
'electromagnetic mass' and no 'material mass'. 

In this connection some known static charged fluid spheres of Tolman- VI type 
solutions have been reexamined and the gravitational masses are shown to be of 
electromagnetic origin by Ray and Das (2007a). They have considered a more 
general form of the gravitational mass as follows: 

2(l + 2n-n2)a ^^"^^^ 

where for physical viability the values on n to be assigned are < n < 2. 

For the specific choice K — l/-\/2 and n = 1 of these parameters, the ansatz 
expressed in equation (1.20) reduces to q{a)/a — l/-\/2, where a is the radius of 
the sphere. It is interesting to note that for this charge-radius ratio all the perfect 
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fluid equations of state reduce to tlie form p + p = wliicli is known as tlie 'pure 
charge condition' (Gautreau 1985) and also imperfect-fluid equation of state in the 
literature for the matter distribution under consideration is in tension and hence the 
matter is named as a 'false vacuum' or 'degenerate vacuum' or 'p-vacuum' (Davies 
1984; Blome and Priester 1984; Hogan 1984; Kaiser and Stebbins 1984). 

Ray and Das (2007b) have again considered the Einstein-Maxwell space-time in 
connection with some of the astrophysical solutions previously obtained by Tolman 
(1939) and Bayin (1978). The effect of charge inclusion in these solutions has been 
investigated thoroughly and the nature of fluid pressure and mass density throughout 
the sphere have also been discussed. Mass-radius and mass-charge relations have 
been found out for various cases of the charged matter distribution. Two cases are 
obtained where perfect fluid with positive pressures gives rise to electromagnetic 
mass models such that gravitational mass is of purely electromagnetic origin. The 
stability conditions have been investigated for all these Tolman-Bayin type static 
charged perfect fluid solutions in connection with the stellar conflgurations. 

1.2.3 Lorentz's Electromagnetic Mass: a Clue for Unifica- 
tion? 

Ray (2007) in his review of the electromagnetic mass model by Lorentz has described 
the philosophical perspectives and given a historical account of this idea, especially, 
in the light of Einstein's Special Relativistic formula E — mc^. It is known that, 
at distances below 10~^^ m, the strong, weak and electromagnetic interactions are 
"different facets of one universal interaction" (Georgy, Quinn and Weinberg 1974, 
Wilczek 1998). This is already confirmed by (i) the theories of the unification of 
electricity and magnetism by Maxwell, (ii) that of earth's gravity and universal 
gravitation by Newton and (iii) "...the unified weak and electromagnetic interac- 
tion between elementary particles..." by S. Glashow, A. Salam and S. Weinberg for 
which Nobel Prize was awarded to them in 1979 . Therefore, as regards unification 
scheme, Ray (2004) has argued that though there has been much progress towards 
a uniflcation of all the other forces - strong, electromagnetic and weak - in the 
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Grand Unified Theory (GUT), gravity has not been included in the scheme. In this 
context it has also been mentioned that there are some problems with gravity: (i) 
the strength of the gravitational interaction is enormously weaker than any other 
force (the hierarchy problem) and (ii) General Theory of Relativity does not consider 
gravity a force, rather a kind of field for which a body rolls down along the space- 
time curvature (the field theoretical problem) . As a probable alternative solution to 
this problem Ray (2004) has put forward Lorentz's conjecture of 'electromagnetic 
mass' and suggested that this may be a competent candidate of the long desired 
unification. 

1.3 Motivation and Discussion of Our Investiga- 
tion 

Our motivation to work on relativistic electromagnetic mass models is based on a 
slightly different approach. We have studied the role of the cosmological constant in 
constructing the electromagnetic mass model. This is the central part of our work. 
It has been observed that Tiwari et al. (1984) and other authors constructed elec- 
tromagnetic mass models without considering any A term. So, we have considered 
here A in the Einstein field equations by assuming it to be a scalar rather than a 
constant and hence re-examined the work of Ray and Ray (1993) and Tiwari and 
Ray (1996). We obtained a class of exact solutions for the Einstein-Maxwell field 
equations by assuming the cosmological constant to be a space variable scalar, i.e., 
A = A(r). The source considered in the chapter II is static, spherically symmetric 
and anisotropic charged fiuid type. The solutions obtained are matched continu- 
ously to the exterior Reissner-Nordstrom solution and each of the four solutions 
represents an electromagnetic mass model. 

In chapter III we have examined whether even without employing the vacuum fluid 
equation of state p + p = a stable model with electromagnetic mass can be con- 
structed. Here we have considered a charged anisotropic static spherically symmetric 
fiuid source of finite radius. The field equations thus obtained under certain mathe- 
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matical assumptions yield a set of solutions which are shown to be electromagnetic 
in origin. Electromagnetic mass models have been studied by several authors under 
the special assumption p + p — 0. Here we have shown that even for p + p 7^ 
electromagnetic mass model can be constructed. This is one of the motivations of 
the present investigation. However, the same question whether there exists any elec- 
tromagnetic mass models where this condition p+p 7^ is violated was addressed by 
Tiwari et al. (1991) and obtained electromagnetic mass model in the isotropic and 
axially symmetric matter distribution or charged dust case only, whereas in chapter 
IV we have searched a solution by employing a relation between the radial and tan- 
gential pressures as p± = p.,. + ax^r^ . Our aim is to see if there is any effect of space 
dependent A on the energy density of classical electron. Here we have considered 
an extended static spherically symmetric distribution of an elementary particle like 
electron having the radius of the order of 10~^^ cm. It is already suggested that 
there might be some negative energy density regions within the particle in the gen- 
eral theory of relativity (Bonnor and Cooperstock 1989, Herrera and Varela 1994). 
It is, therefore, argued in the present investigation that such a negative energy den- 
sity also can be obtained with a better physical interpretation in the framework of 
Einstein-Cartan theory. 

In many theories higher dimensions play an important role, specially in superstring 
theory which demands more than usual four dimensional space time. This is also 
true in studying the models regarding unification of gravitational force with other 
fundamental forces in nature. So long electromagnetic mass model has been stud- 
ied extensively in the four dimensional space-times of the General Relativity. Here, 
in chapter VI, we have presented electromagnetic mass model in the space-time of 
higher dimensional theory of general relativity. Under this motivation we have con- 
sidered here a static spherically symmetric charged dust distribution corresponding 
to higher dimensional theory of general relativity. 

In chapter VI we have studied static spherically symmetric anisotropic source for 
the Einstein- Maxwell space-times assuming the erstwhile cosmological constant A 
dependent on the spatial coordinate, viz., A = A(r). It is shown that the solutions 
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thus obtained are of electromagnetic in origin in the sense that all the physical pa- 
rameters including the gravitational mass originate from the electromagnetic field 
alone. It is also shown that the generally used pure charge condition, viz., p-\-pr — 
is not always required for constructing electromagnetic mass models. 
The concluding chapter VII offers a general discussion on the whole work along with 
the future scope of the field of electromagnetic mass models in General Relativity. 
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Chapter 2 



Relativistic Electromagnetic Mass 
Models with Cosmological 
Variable A in Spherically 
Symmetric Anisotropic Source 



"Whether one or the other of these methods will lead to 
the anticipated "world law" must be left to future research. 
- Max Born (1962) 



2.1 Introduction 

A very important problem in cosmology is that of the cosmological constant the 
present value of which is infinitesimally small (A < 10~^^ cm~^). However, it is 
believed that the smallness of the value of A at the present epoch is because of 
the Universe being so very old (Beesham 1993). This suggests that the A can not 
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be a constant. It will rather be a variable, dependent on coordinates - either on 
space or on time or on both (Sakharov 1968; Gunn and Tinslay 1975; Lau 1985; 
Bertolami 1986a,b; Ozer and Taha 1986; Renter and Wetterich 1987; Preese et al. 
1987; Peebles and Ratra 1988; Wampler and Burke 1988; Ratra and Peebles 1988; 
Weinberg 1989; Herman et al. 1989; Chen and Wu 1990; Berman and Som 1990; 
Abdel- Rahman 1990; Berman 1990a,b; Berman 1991a,b; Sistero 1991; Kalligas et 
al. 1992; Carvalho et al. 1992; Ng 1992; Beesham 1993 and Tiwari and Ray 1996). 
Now, once we assume A to be a scalar variable, it acquires altogether a different 
status in Einstein's field equations and its infiuence need not be limited only to 
cosmology. The solutions of Einstein's field equations with variable A will have a 
wider range and the roll of scalar A in astrophysical problems will be of as much 
significance as in cosmology. 

It is this aspect that motivated us to reexamine the work of Ray and Ray (1993) 
and Tiwari and Ray (1996) with the generalization of anisotropic and charged source 
respectively. One can realize from the present investigations how the variable A gen- 
erates different types of solutions which are physically interesting as they provide a 
special class of solutions known as electromagnetic mass models (EMMM). 
In section 2.2, the Einstein- Maxwell field equations with variable A are derived. So- 
lutions corresponding to different cases for anisotropic system are obtained in section 
2.3. All the solutions obtained are matched with the exterior Reissner-Nordstrom 
(RN) solution on the boundary of the charged sphere. Pinally, some concluding 
remarks are made in section 2.4. 

2.2 Field Equations 

The Einstein-Maxwell field equations for the spherically symmetric metric 

ds^ = e'^^'Ut'' - e^^'Ur^ - r^(de^ + sin'^edcj)^) (2.1) 
corresponding to charged anisotropic fiuid distribution are given by 

e-\X'/r - l/r") + l/r^ = Sirp + E'' + A, (2.2) 
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e-^{y' jr + l/r^) - l/r^ = 87rp^ - E"^ - (2.3) 

e-^[i/72 + i/'74 - i/'A74 + (i/' - A') /2r] = Sttp^ + - A (2.4) 

and 

(r^^)' = 47rrVe^/l (2.5) 

The equation (2.5) can equivalently be expressed in the form, 

Eir) = \ r Anr^ae'^/^dr = (2.6) 

where q{r) is total charge of the sphere under consideration. 
Also, the conservation equation is given by 

^{pr - A/Stt) + {p + Pry/2 = g^|:(g') + 2(p± -Pr)/r. (2.7) 

Here, p, p^, p_L, E, a and q arc respectively the matter-energy density, radial 
and tangential pressures, electric field strength, electric charge density and electric 
charge. The prime denotes derivative with respect to radial coordinate r only. 
Equations (2.2) and (2.3) yield 

e-^(t/' + A') = 87rr(p + p,). (2.8) 

Again, equation (2.2) can be expressed in the general form as 

e-^ = l-2M(r)/r, (2.9) 

where 

Mir) = 4:71 Hp+iE^ + A)/87iydr (2.10) 

^0 

is the active gravitational mass of a charged spherical body which is dependent on 
the cosmological parameter A = A(r). 
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2.3 Solutions 



A number of solutions can be obtained depending on different suitable conditions 
on equation (2.7). Here we assume the relation (j'oofi'ii = —1, between the metric 
potentials of metric (2.1), which, by virtue of equation (2.8), is equivalent to the 
equation of state ^ 

p+Pr = Q. (2.11) 

The equations (2.2) - (2.5) being underdetermined, we further assume the following 
conditions 

o-e^/^ = (To (2.12) 

and 

P± = npr, {n 7^ 1), (2.13) 

where (Tq is a constant (which from (2.6) can be interpreted as the volume density of 
the charge a being constant) and n is the measure of anisotropy of the fluid system. 
Equation (2.6), with equation (2.12), provides the electric field and charge as 

E = q/r'^ = 47r(7or/3. (2.14) 

Using equations (2.11), (2.13) and (2.14), in equation (2.7), we get 

-^{Pr-A/8Tr)-2{n-l)pr/r = 2Ar, A = 2na^/3, (2.15) 
ar 

which is a linear differential equation of first order. 

Since the equation (2.15) involves two dependent variables, Pr and A, to solve 
this equation, we consider the following four simple cases. 



■^In terms of energy- momentum tensor this can be expressed as T°o = T^i. 



27 



2.3.1 A = Ao — 87rpr, (Aq =constant) 

The solutions in this case are obtained as 



e'^ = e^^ = 1 - 2M(r) /r, (2.16) 
p=-Pr = -p±/n = (A - Ao)/87r = ^a-^^-'V^K'^ - r"-=']/(n - 3) (2.17) 

and 

M(r) = %w ^^J (ri + 2)(n + 3)a"-^ - SOr^-^] + Aor^G, (2.18) 

where a is the radius of the sphere. 

Some general features of these solutions are as follows: 
(1) As we want, customarily, p > (and hence Pr < 0), we must have, from (2.17), 
n > 3. However, we can choose n < 3 (and certainly n ^ 1). In that case also 
p becomes positive. This result, viz., the positivity of matter-energy density is 
obvious as the electron radius for the present model is 10~^^ cm, which is much 
larger than the experimental upper hmit 10~^^ cm (Quigg 1983). Within this limit 
the charge distribution of matter must contain some negative rest mass (Bonnor 
and Cooperstock 1989; Herrera and Varela 1994). This is the reason why we cannot 
consider p < and hence pj. > in the present model. 

(2) Similarly, we can observe that the effective gravitational mass (which we get 
after matching of the interior solution to the exterior RN solution on the boundary) , 

m = M(a) + g2(a)/2a - AoaVS = SnA{n + 3)aV5(n + 2), (2.19) 

is positive for both the choices, n > and n < 0. In this respect, the Tolman- 
Whittaker mass, 

rriTw = / {T% - T'^o)yf^dV, (a = 1, 2, 3 and AD) 
87r/lr7~("~^)r^ 

[2(n + 2)(n + 3)a"-3-15(n+l)r"-3] -AorV3, (2.20) 



15(n-3)(n + 2) 
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can also be examined. In general, this is negative and also equal to modified Tolman- 
Whittaker mass (Devitt and Florides 1989), 

moF = e-^''+^~^'^mTw, (2.21) 

as + A = 0, by virtue of the condition gooQn = — 1 in the present paper. 
(3) Pressure being negative the model is under tension. This repulsive nature of 
pressure is associated with the assumption (2.11), where matter-energy density is 
positive. This negativity of the pressure corresponds to a repulsive gravitational 
force (Ipser and Sikivie 1983; Lopez 1988). 

(4) The cosmological parameter A, which is assumed to vary spatially, can be 
shown to represent a parabola having the equation of the form A = 87r>l[(a/2)^ — 
(r — a/2)^] + Aq for a particular case n — 2. The value of A increases from to 
a/2 and then decreases from a/2 to a and hence it is maximum at a/2 . The vertex 
of the parabola is at r = a/2 whereas the values of A at r = and at r = a are 
Aq, the erstwhile cosmological constant. The same result can also be obtained from 
equation (2.17) as at the boundary of the sphere r — a, — p± — Q { and hence 
A = Ao). 

(5) The solution set provides electromagnetic mass model (EMMM) (Feynman 
et al. 1964; Tiwari et al. 1984, 1986, 1991; Gautreau 1985; Gr0n 1985, 1986a, 1986b; 
Ponce dc Leon 1987a, 1987b, 1988; Tiwari and Ray 1991a, 1991b, 1997; Ray et al. 
1993; Ray and Ray 1993). This means that the mass of the charged particle such 
as an electron originates from the electromagnetic field alone (for a brief historical 
background, see Tiwari et al. 1986). 

(6) The present model corresponds to Ray and Ray (1993) for n — 1, under 
the assumption — — A/Stt. It can be observed that the other simple possibility, 
Pr — A/Stt, does not exist for this case (equation (23) of Ray and Ray (1993)). 
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2.3.2 A = Ao + STrpr 



In this case we have the following set of solutions: 




l-2M(r)/r, 



(2.22) 



■p±/n 



(A- Ao)/87r = ArV(n- 1) 



(2.23) 



and 



M(r) = 47rArVl5 + AorV6. 



(2.24) 



Here some simple observations are as follows: 

(1) In this case also the electron radius being ~ 10~^^ cm the matter-energy density 
should be positive (Bonnor and Cooperstock 1989; Herrera and Varela 1994). This 
positivity condition requires that n must be greater than unity. 

(2) The effective gravitational mass, 



is always positive whereas the Tolman-Whittaker mass which is also equal to the 
modified Tolman-Whittaker mass, i.e.. 



is always negative in the region < r < a. The gravitational mass in this case is 
independent of anisotropic factor n. 

(3) The pressures pj. and pj_ are repulsive for n > 1 (as in the previous case). 

(4) The equation (2.23) for n = 2 can be written in the form A = —SirAr'^ + Aq. 
This yields a half-parabola whose vertex is at r = and the parabola lies in the 
fourth-quadrant of the coordinate systems (r. A). 

(5) The effective gravitational mass as obtained in (2.25) is of electromagnetic ori- 
gin. 

(6) The matter-energy density p as well as the pressures Pr and p± are all zero at 
the centre of the spherical distribution and increase radially being maximum at the 



/5, 




rriTw — 'iTiDF — —IQnAr^ /lb — Aor'^/3, 



(2.26) 
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boundary. This situation is somewhat unphysical though not at all unavailable in 
the literature (Som and Bedran 1981). 



2.3.3 A = Ao-87r/^(/r 

The solution set for this case is given by 

= e-^ = 1 - 2M(r) /r, (2.27) 
p = -p^ = -pjn = 2yla-(2"-5VV(2n - 5)[a2"-^ - r^''-% (2.28) 

A = (2T- 3)(2n-5) ^^^"^ ~ " ^""'""'^ " 87rAa7(2n - 3) + Ao (2.29) 

and 

^^<'-' = 15„(2n-3)(2n-5) '"'" + " " ^^<" " I"-""' 

-47rAaVV[3(2n - 3)] + K^r^/Q. (2.30) 

Some general features of the above set of solution arc as follows: 

(1) The matter-energy density is positive and pressures are negative for n > 5/2. 

(2) The effective gravitational mass, 

m = 87rA(3n + l)a^ /Ihn, (2.31) 

is positive for n > 1. On the other hand, the Tolman-Whittaker mass and the 
modified Tolman-Whittaker mass, being equal, are given by 

327i-y4a-(2"-5)^5 

mr^ = m.. = - i5,(2n-3)(2n-5) ["^" + '^^'^ " 

-15(n - 1) (2n - l)r2"-^ + 87r^aVV[3(2n - 3)] - AorVS. (2.32) 

Depending on the different values of n these masses may be negative or positive. 

(3) The matter-energy density and the pressures, as usual, are zero at the centre 
r = as well as at the boundary r — a. Thus the maximum value must be in the 
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region < r < a. This can be confirmed from equation (2.28) whicli, for the value 
n — 2, represents a parabola of the form A = 2A[(a/2)^ — (r — a/2)^], the vertex 
being at r = a/2. 

(4) The value of A at the centre r = is [Aq — 87rAa^/(2n — 3)]. It acquires 
maximum value Aq at the boundary r — a. 

(5) The solution set represents EMMM. 



2.3.4 A = Ao + /^(ir 

The solutions in this case are given by 

e^ = e"^ = l-2M(r)/r, (2.33) 

P=-Pr = -vUn = '-^ --\a^-^ - r^-\ (2.34) 

^ = - (2:-l)(2,.-3) 'P"-^"""'-^-""i 

+87rAaV(2n - 1) + Aq (2.35) 



and 



^ - 15(. + l)(2.-l)(2.-3) ^-(- + - - ^^^^ - 

Wr3 Aor3 
3(2n-l) 6 ^ ' 

Here, the observations are as follows: 

(1) The matter-energy density is positive and pressures are negative for n > 3/2. 

(2) The effective gravitational mass, 

m = 87rA(3n + 5)aVl5(n + 1), (2.37) 

for the condition n > 1 is always positive, whereas the Tolman-Whittaker mass, 

87rAa-(2"-3)r5 

m.^ ^mn.^- ^^^^ ^ ^^^^^ _ ^^^^^ _ [4n(n + l)(2n - l)a " 

-15(n - l)(2n + l)r2"-3] - 87rAaVV[3(2n - 1)] - AorV3 (2.38) 



32 



may have positive or negative value depending on the choice of n. 

(3) The values related to p and p are zero both at r = and r — a. 

(4) The effective gravitational mass as well as the other physical variables, including 
A, are of purely electromagnetic origin. 



2.4 Conclusions 

(1) As mentioned in the introduction, the present work considers A, the erstwhile 
cosmological constant, to be a variable dependent on space coordinates. The con- 
tribution of this variable A can be seen in the calculations given in the previous 
sections. It can be seen that A is related to pressure and matter-energy density, and 
therefore contributes to effective gravitational mass of the astrophysical system. 

(2) The present EMMMs have been obtained under the condition p + pr — (equa- 
tion (2.11)). This problem thus requires further investigation to see whether such 
models can be obtained even for the condition p + Pr 7^ 0. 



The contents of this chapter published in Indian Journal of Pure and Applied 
Mathematics (2000) 31 1017. 
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Chapter 3 



Classical Electron Model with 
Negative Energy Density in 
Einstein- Cart an Theory of 
Gravitation 



"// you can look into the seeds of time, and say 
which grain will grow and which will not. . . " 
- Shakespear {Macbeth) 



3.1 Introduction 

Recently, Cooperstock and Rosen (1989), Bonnor and Cooperstock (1989), and Her- 
rera and Varela (1994) have shown that within the experimentally obtained upper 
hmit of the size of the electron ( ~ 10~^^ cm) (Quigg 1983), when it is modeled 
as a charged sphere obeying Einstein-Maxwell theory, must contain some negative 
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gravitational mass density regions within the particle . According to Cooperstock, 
Rosen and Bonnor (CRB) (1989), the rest mass or active gravitational mass within 
this sphere, by virtue of the relation 

M = ™-|1, (3.1) 

is negative and about 10~^^ cm (when the inertial mass or effective gravitational 
mass, charge and radius, respectively, of the electron, are m — 6.76 x 10~^^ cm, 
q — 1.38 X 10~^^ cm, and a — 10"^^ cm in relativistic units). Further, Herrera 
and Varela (HV) (1994) have shown, in one of the cases of their paper, that the 
matter-energy density 

P=(«g' + ^7rao')(a'-r2), (3.2) 
for the constant a = —4.77 x lO^^cm"^ (when radius a ~ 10^^^ cm) is also negative, 
(To being the constant charge density at the centre of the spherical distribution. 
These models, however, lack spin and magnetic moment and hence do not possess 
the actual physical characteristics required for an electron. 

As an alternative way both the groups suggest the stationary Kerr-Newman (KN) 
metric (Newman et al. 1965) related to the solution of Einstein-Maxwell equations 
to be more appropriate than those described earlier. However, in this context it is 
also to be mentioned here that the KN metric cannot be valid for distance scales 
of the radius of a subatomic particle (Mann and Morris 1993; Herrera and Varela 
1994). 

We, therefore, feel that the problem can be tackled in the framework of Einstein- 
Cartan (EC) theory, where torsion and spin are inherently present in the formulation 
of the theory itself. 



3.2 An Overview: The Negative Density Models 

Before going into the Einstein-Cartan theory let us have a bird's-eye view of the 
negative matter-energy density models which we have mentioned in the introduc- 
tion. 
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3.2.1 The Cooperstock-Rosen-Bonnor (CRB) Model 

Cooperstock and Rosen (1989) and Bonnor and Cooperstock (1989) in their papers 
have shown that any spherically symmetric distribution of charged fluid, irrespective 
of its equation of state, whose total mass, radius and charge correspond to the 
observed values of the electron, must have a negative energy distribution (at least 
for some values of the radial coordinate). Considering a static spherically symmetric 
charge distribution with the line-element 

ds-" = e^^^'Uf - e^^^'^dr^ - r^dO^ + sin^edcl)^) (3.3) 

they have argued that when the Einstein- Maxwell equation 

R% - b%R = 87r(r%^"^^ + r%^''"^) (3.4) 



is written in the form 



1 - - r(87rp + e-^'^^^^E^ydr (3.5) 
r Jo 



and hence is equated with the Reissner- Nordstrom exterior metric on the boundary 
r — a, which as usual gives 

1 _ ^ + 4 = 1 - 1 r{87rp + e-^'^+'^^Ey^dr. (3.6) 
a a Jo 



Then for the previous specifications of mass, charge and radius of the electron it can 
be shown that 

4 - — ~ 2 X 10-36 > 0. (3.7) 

So, the left hand side of the above equation (3.6) must be greater than unity and 
hence on the right hand side p < for some values of r implying that the electron 
must contain some negative rest mass density though the net mass is as usual a 
positive quantity. 
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3.2.2 The Herrera-Varela (HV) Model 

Following the CRB model (1989) Herrera and Varela (1994) have discussed the 
fact that the electron, when modeled as a relativistic spherically symmetric charged 
distribution of matter, must contain some negative rest mass if its radius is not 
larger than ~ 10^^^ cm. In this regard they have analyzed some extended electron 
models and have shown that negative energy density distributions result from the 
requirement that the total mass of these models remains constant in the limit of a 
point particle. Among all these extended electron models the model of Tiwari et 
al. (1984) demands special attention to us which will be seen very much relevant to 
our present work. Herrera and Varela (1994) generahze this model of Tiwari et al. 
(1984) by introducing a condition of anisotropy in the form 

P± — Pr — Oi(fT^ (3.8) 

where a is a constant. 

Thus the solution obtained by Herrera and Varela (1994) is as follows: 

45 15 
Pr = -{aq^ + '^Ti<Jo^){a^ - r'), (3.10) 
= aq'r^ - {aq^ + ^7r(7o')(a2 - r'), (3.11) 

O 

64 8 

m = — TT^o-Q^a^ H Trag^a^, (3.12) 

45 15 

P= («g' + ^W)(a'-r2) (3.13) 

and 

q = ^-KGoa^ (3.14) 

The value of a can be obtained from the equation (3.12) as a = —4.77 x 10^^ cm^^ 
and hence the energy density, as given by the equation (3.13) is negative for the 
radius of the electron a = 10~^^ cm. 

Now, from the equation (3.12) it can be seen that the effective gravitational mass. 
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m, is of purely electromagnetic origin and corresponds to the TRK model (1984) 
with a — case. This type of models where mass, including all the other physical 
parameters, originate from the electromagnetic field alone are known as the electro- 
magnetic mass models [EMMM] in the literature [Feynman et al. 1964] and have 
been investigated by several authors (Florides 1962, 1983; Cooperstock and de la 
Cruz 1978; Tiwari et al. 1984, 1986, 1991, 2000; Gautreau 1985; Gr0n 1985, 1986a,b; 
Ponce de Leon 1987a,b, 1988; Tiwari and Ray 1991a,b, 1997; Ray et al. 1993; Ray 
and Ray 1993). In the present paper we shall construct such a model within the 
framework of Einstein-Cartan theory with negative matter-energy density for some 
values of the radial coordinate. 



3.3 The Field Equations of Einstein-Cartan The- 
ory 

The EC field equations are given by 

R'j-l3'jR-->^t'j (3.15) 

and 

Q'jk - S'jQ\k - S'kQ'ji = ->^s'j,, (3.16) 

where fj is the canonical energy- momentum tensor (asymmetric), Q^jj^ is the torsion 
tensor and S^jk is the spin tensor (with k — —Sir, G and c being unity in relativistic 
units) . 

The asymmetric energy-momentum tensor here is given by 

fj^T'j + lg''Vm{S"'jk), (3.17) 

being covariant derivative with respect to the torsionless, symmetric Levi-Civita 
connection and the symmetric energy-momentum tensor T'^ will consist of two 
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parts, viz., matter and electromagnetic tensors and which, respectively, are 

T-"'^ = {p + p)u'u,-pg'^ (3.18) 

and 

^i[em) ^ 1_^_F^^F''' + -^5',FuiF^'), (3.19) 

where p is the matter-energy density, p is the fluid pressure, is the velocity four- 
vector (with u^Ui — 1) and Fij is the electromagnetic field tensor. 
The conservation equations for the EC theory can be given through the Bianchi 
identities as 

Vk[{p + pV-g'''u'Vm{u'^Su)]^u^V^p (3.20) 

and 

[{p+p)u^-g^'u'Vm{u^S'')]Vi,u^ = -Vi{u%)+v!'SimR'^k-]^u^Si^R'^jk (3.21) 

Now, electromagnetic fields not being coupled with torsion (Novello 1976; Raychaud- 
huri 1979) the Maxwell equations as usual take the form 

VjF'^ = f (3.22) 

and 



{J'^),i = 0. (3.23) 

The electromagnetic field tensor, F^, in the above equation (3.22) is related to the 
electromagnetic potentials as = Ajj — Aj^i which is equivalent to ^.j = 0, Ai 
being the electrostatic potentials. Here and in what follows a comma denotes the 
partial derivative with respect to the coordinate indices involving the index. 

Again, the spin tensor and the intrinsic angular momentum density tensor are 
related in the form 

S',k = u'Sjk, (3.24) 

with 

u'Sik = 0. (3.25) 



39 



Now, assuming that the spins of the individual charged particles composing the fluid 
distribution are all aligned in the radial directions (Prasanna 1975; Raychaudhuri 
1979; Tiwari and Ray 1997) and the matter is at rest with respect to the observer, 
the non-vanishing components of the spin tensor can be obtained, from equations 
(3.24) and (3.25), as 

5°23 = -^°32 = s{g,,)-'l\ (3.26) 
whereas, from equation (3.16), we have the torsion tensor as 

Q\z = -Q°32 = -^s{go^)-'/\ (3.27) 

s = S23 being the only non vanishing component of the intrinsic angular momentum 
density tensor. Here, we have followed the convention {t, r, 9, (j)) — (0, 1, 2, 3). 

The Einstein-Cartan-Maxwell equations with source can be written as (Tiwari 
and Ray 1997) 



+ + (3.28) 



e--(^ + l)-l = 8.p~-£;^ (3.29) 



\v" + - I/' A' + ^ ^] = SttpI + ^2 (3.30) 

r 

and 

{r'^E)' = 47rrVe^, (3.31) 

where p, Rl and E are the effective matter-energy density, effective pressures 
(radial and tangential) and electric field respectively, and are defined as 

p = p - 27rs^ (3.32) 

Pr — Pr ~ 27rs^, (3.33) 

p±^p±- 2tts^ (3.34) 

and 

E = -exp[-{iy + A)](/.' = ^, (3.35) 

Pr, P±, s, (f) and q being the ordinary matter-energy density, ordinary pres- 
sures (radial and tangential) , spin density, electrostatic potential and electric charge 
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respectively. Here, a represents the electric charge density and prime denotes dif- 
ferentiation with respect to the radial coordinate r. 

Now, the conservation equations (3.20) and (3.21) with the help of equations (3.22) 
and (3.23) reduce to 

— = -(p + p.). + ^-(g ) + . (3.36) 

This is the key equation which is to be solved for constructing EMMM. 



3.4 The Solutions 

Addition of (3.28) and (3.29), under the assumption g'ooS'ii = — 1 ( or equivalently, in 
terms of energy- momentum tensors T°o — provides the pure charge condition 

p + p,. = 0, (3.37) 

where, in general, p is assumed to be positive and hence p,. is negative. However, 
as is evident from equation (3.32), p, being the effective energy-density, can even be 
negative due to the positive second term related to the spin on the right hand side. 
Thus, the possibility of equation (3.33) being satisfied with being positive is not 
ruled out. 

To make (3.31) and (3.36) solvable, we further assume that 

cre^ = (7o (3.38) 

and 

P± — = P_L — = ac^r^ (3.39) 

following TRK Model (Tiwari et al. 1984) and HV model (Herrera and Varela 1994) 
respectively, where ctq and a are two constants as mentioned earlier. 
By substituting (3.37) - (3.39) in (3.31) and (3.36), we get 

^ = 4 = ^^'^or, (3.40) 
3 
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and 



Pr = 2ns^ - {aq^ + '^W){a^ - r^), (3.41) 

2 

p± = 27rs^ - aq^ia'^ - 2r'^) - -7rc7o^(a^ - r^) (3.42) 



p = 27rs^ + (ag^ + ^TTCTo^) (a^ - r^) . (3.43) 



The active gravitational mass 

Mir) = Att Hp - 2ns^ + ^^dr (3.44) 

Jo OTT 

takes the form, by virtue of (3.40) and (3.43), as 

M(r) = 4-^^<^o^r2[87rQ;a^(5a2 - 3r^) + 3(5a^ - 2r% (3.45) 
135 

Thus, the metric potentials A and u are given by 

= = 1 _ ?MM, (3.46) 



r 



whereas, the effective gravitational mass mentioned in (3.1), can be obtained as 

64 2 
m = ^7rVo2a5(l + ^Traa^), (3.47) 
45 3 

which corresponds to the second case (B) of HV model (1994) and is of purely 
electromagnetic origin. This corresponds to the TRK model (1984) with a = case. 
It can be noted here that unlike the matter-energy density the effective gravitational 
mass is independent of spin. 

In this context it is to be mentioned here that the junction conditions in the EC 
theory are different from that of general theory of relativity and indeed read like 
this (Arkuszewski et al. 1975) 

mu'l^^O (3.48) 

and 

p |_= 27rG{niS') |_ (3.49) 

where 5"* is the spin density pseudo- vector. Here condition (3.48) is the same as in 
classical relativistic hydrodynamics and has already been incorporated by matching 
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the interior solution with the exterior Reissner-Nordstrom field at the boundary of 
the spherical distribution. The condition (3.49), however, is the additional condition 
to be satisfied in EC theory. In the present case, it is only the effective pressure 
(radial) that vanishes on the boundary and not the ordinary radial pressure which, 
by virtue of equation (3.41), equals 27rs^. The spin is ahgned in the radial direction 
and hence the spin density pseudo- vector is hypersurface orthogonal. Thus, the 
boundary condition (3.49) will become 

Pr\r=a-G = 27rs\^„_o, (3.50) 

which, depending on whether s is a constant or function of coordinates, will auto- 
matically be satisfied. 

In this connection it is to be mentioned here that the spin density, s, in the final 
solutions (3.41) - (3.43) remains arbitrary (function of r). An explicit functional 
form of this spin density can also be obtained by assuming some additional phys- 
ically viable possibility, such as the one used by Prasanna (1975) by splitting the 
conservation equation into two parts, the second part relating to conservation of 
spin only, giving the functional form of spin density as s = SQe~'' (where Sq is the 
value of s at r = i.e. the central spin density). This can, using equations (3.38) 
and (3.46) and the condition Qoogu — 1, (that is, iz + X — O) equivalently be written 
as S(T = sqCo = constant. The functional form of spin density is, however, not rele- 
vant in our discussion as our problem is concerned with the properties related to the 
'electron', an elementary particle whose radius is of the order of 10~^^ cm. Indeed, 
as the spin function is arbitrary, there is no loss in generality, even if we assume it 
to be almost a constant (that is, the quantized value of the spin of the electron). 

3.5 The Negative Energy Density Model 

Let us have a closer observation of the results of the previous section 3.4. The 
equation (3.43) related to matter-energy density has the spin density part in the first 
term where spin density is defined as s = SS/Aira^, where S is the spin of electron 
the quantized value of which is S — h/2. Then substituting the standard values for 
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different parameters in the relativistic units, as mentioned in the introductory part, 
the numerical value for the matter-energy density (3.43) can be shown as 

p = 6.14 X 10-^^ - 6.81 X 10^^(10-^^ - r^). (3.51) 

The first term related to spin, being of the order of 10~^^, is too small compared to 
the last term and hence the spin contribution is negligible. Now, the equation (3.51) 
indicates that the central density at r = is negative and its magnitude is about 
10~^. On the other hand, the total density at the boundary, r = a, is positive as 
usual with the numerical value about 10~^^. This change in the sign of the energy 
density is because of the presence of the spin term in equation (3.51) which, indeed, 
is the contribution of the EC theory. In the absence of spin, however, we could have 
negative and zero densities at the centre and boundary of the electron respectively. 
This change in the sign again indicates that the central negative value gradually 
increases along the radius and somewhere, in the region < Tc < a, it becomes zero, 
where Tc is the critical radius. Obviously, the amount of negative energy density is 
less than its positive counterpart the balance of which ultimately provides the net 
density as the positive one. 

It is already mentioned that, in general, for any spherical fluid distribution the 
density on the surface should be zero where we are getting some non-zero value for 
it. This finite value of density is solely coming from the spin contributed part 27rs^. 
Thus for the vanishing spin the situation corresponds to the general behaviour (vide 
equation (17) of Herrera and Varela 1994). In this context it is also to be noted here 
that up to the critical radius behaviour of our model is similar to that of Herrera 
and Varela (1994). Beyond this cut off radius the energy density is regulated by spin 
which makes the overall density of the model positive. This particular aspect lack in 
the model of Herrera and Varela (1994) where the total energy density is a negative 
quantity. This increase of matter-energy density due to spin density can probably 
be accounted for the kinetic energy through the angular motion of the electron here. 
Similar kind of examination is also possible for the pressures, radial and tangential, 
both. The radial pressure, in this case of equation (3.41), takes the following value 

Pr = 6.14 X 10-^^ + 6.81 X 10^^(10-=^^ - r^). (3.52) 
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However, pressure is throughout positive here from the centre to boundary. These 
results, i.e. negative energy density and positive pressure, are in accordance with 
the pure charge condition (3.37) which reads as p = —pr + 47rs^ . 

3.6 Conclusions 

(i) The possible origin of the intriguing negative matter-energy density in the work 
of Cooperstock and Rosen (1989), Bonnor and Cooperstock (1989), Herrera and 
Varela (1994) and present paper may be due to the finiteness of the total mass of 
the Reissner-Nordstrom solution (Visser 1989). Since the electrostatic energy of a 
point charge is infinite, the only way to produce a finite total mass is the presence 
of an infinite amount of negative energy at the center of symmetry. According to 
Bonnor and Cooperstock (1989) the negativity of the energy density and hence the 
active gravitational mass is consistent with the phenomenon, known as the Reissner- 
Nordstrom repulsion (de la Cruz and Israel 1967; Cohen and Gautreau 1979; Tiwari 
et al. 1984; Cooperstock and Rosen 1989). In this regards Bonnor and Cooper- 
stock (1989) also have discussed about the singularity theorems of general relativity 
(Hawking and Ellis 1973). They have shown that the negative regions are hable 
to exist over distances of order 10~^^ cm and as the proof of the singularity theo- 
rems depends on the manifold structure of space-time vahd down to lengths of order 
10~^^ cm so might break down below this. On the other hand, in the context of 
Einstein-Cartan theory of gravitation the idea of negative mass is not a new one as 
stated by de Sabbata and Sivaram (1994): "...torsion provides a natural framework 
for the description the negative mass under extreme conditions of such as in the 
early imiverse, when a transition from positive to negative mass can take place." 

(ii) We have considered in the present chapter an extended static spherically sym- 
metric distribution of an elementary particle like electron having the radius of the 
order of 10~^^ cm, and even if for the finite size of the physical system the spin in 
Einstein-Cartan theory can be related to orbital rotation (which indeed is not the 
case), for systems of dimensions of subatomic particle the orbital rotation loses its 
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meaning. In this case, the only way is to take the spin to be the 'intrinsic angular 
momentum', that is, the spin of quantum mechanical origin (in our problem since s 
is arbitrary, we can consider its quantized value or an average value) . In this respect 
we would like to quote here from Hehl et al. (1974), "It is crucial to note that spin 
in C/4 theory is canonical spin, that is, the intrinsic spin of elementary particles, not 
the so-called spin of galaxies or planets." 

(iii) Following other authors (Prasanna 1975; Raychaudhuri 1979; Tiwari and Ray 
1997), in the present work the spins of all the individual particles are assumed to 
be oriented along the radial axis of the spherical systems. As to how this alignment 
is brought about is not very much clear. We have discussed here only a few possi- 
ble ways of realizing this situation. According to Raychaudhuri (1979), in general, 
there will be an interaction between the spins of the particles and the magnetic field. 
The overall effect is the alignment of the spins. In this context, Prasanna (1975) 
mentioned that such an alignment may be meaningful either in the case of spherical 
symmetry when magnetic field is present or else one has to consider axially sym- 
metric field. As stated above, our view point is that in the case of physical systems 
of the size of the electron, the radial alignment of the spin is not ruled out. The 
solution obtained supports this view. 

(iv) Though our present approach via Einstein-Cartan theory to inject spin may be 
interesting, we feel even that there should have some room to discuss the relation- 
ship of our work with an alternative means to provide spin and magnetic moment. 
We think this may possible through Dirac-Maxwell theory where spin and magnetic 
moment are naturally incorporated through the Dirac spin (Bohun and Coopcrstock 
1999; Lisi 1995) and would like to pursue this problem in future investigations. 



The contents of this chapter published in International Journal of Modern Physics 
D (2004) 13 555. 
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Chapter 4 



Energy Density in General 
Relativity: a Possible Role for 
Cosmological Constant 



"The most insignificant thing contains some 
little unknown element. We must find it!" 
- Maupassant 



4.1 Introduction 

The structure of electron was, for a long time, an intrigue problem to the researchers. 
Many scientists, like Lorentz (1904) and even Einstein (1919) tried to solve the prob- 
lem to show that the electron mass is a electromagnetic field dependent quantity 
(for a detail account see the references Tiwari, Rao and Kanakamedala (1986) and 
Wilczek (1999)). Later on, under general relativity some models have been con- 
structed by different authors describing extended electron with its mass entirely of 
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electromagnetic origin (Florides 1962; Coopcrstock and dc la Cruz 1978; Tiwari, 
Rao and Kanakamedala 1984; Gautreau 1985). Recently, based on the experimen- 
tal upper limits on the size of the electron as ~ 10~^^ cm (1983) it is argued by 
Cooperstock and Rosen (1989), Bonnor and Cooperstock (1989) and Herrera and 
Varela (1994) that in the framework of general theory of relativity the electron- like 
spherically symmetric charged distribution of matter must contain some negative 
mass density. Being motivated by these results with historical and heuristic values 
we would like to explore a possible role for cosmological constant on the energy 
density of electron when it is modeled as a dependent on the radial coordinate r of 
the charged spherical matter distribution. 

The basic logic for considering variability of so called cosmological constant, which 
was introduced by Einstein in 1917 to obtain a static cosmological model, is related 
to the observational evidence of high redshift Type la supernovae (Perlmutter et 
al. 1998; Riess ct al. 1998) for a small decreasing value of cosmological constant 
{h-present < 10~^^cm~^) at the present epoch. This indicates that instead of a strict 
constant the A could be a function of space and time coordinates. If the role of time- 
dependent A is prominent in the cosmological realm, then space-dependent A has 
an expected effect in the astrophysical context. It is, therefore, argued by Narhkar 
et al. (1991) that the space-dependence of A cannot be ignored in relation to the 
nature of local massive objects like galaxies. Our aim, however, to see if there is any 
effect of space-dependent A on the energy density of the classical electron. This is 
because cosmological constant is thought to be related to the quantum fluctuations 
as evident from the theoretical works by Zel'dovich (1967). Moreover, it is believed 
through indirect evidences that 65 % of the contents of the universe is to be in the 
form of the energy of vacuum (Martins 2002). Thus, the energy density of vacuum 
due to quantum fluctuation might have, in our opinion, some underlying relation to 
the energy density of Lorentz's extended electron (1904) under general relativistic 
treatment. 

In the present chapter IV we have tried to flnd out, through some speciflc case stud- 
ies, that energy density of classical electron is related to the variable cosmological 
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constant and the gravitational mass of the electron is entirely dependent on the 
electromagnetic field alone. 

4.2 The field equations 

To carry out the investigation we have considered the Einstein-Maxwell field equa- 
tions for the case of anisotropic charged fiuid distribution (in relativistic units 
G — c — 1) which are given by 

G'j = R'j - g'jR/2 = -STrfr/""^ + ^V'""^ + rV'"'^], (4.1) 

[{-9f"F% = A7rr{-gf' (4.2) 

and 

F[ij,k] = (4.3) 

where F*-' is the electromagnetic field tensor and J*, current four vector which is 
equivalent to J* = (ru\ a being the charge density and is the four- velocity of the 
matter satisfying the relation Uiu^ ~ 1. 

The matter, electromagnetic and vacuum energy-momentum tensors are, respec- 
tively given by 

r/™^ = (P + PA.)u'uj - p±g'j + (px - PrWVj, (4.4) 
j.i,(em) ^ _^p.^pik _ gi.F^^F^i/A]/AT: (4.5) 

and 

r/^"'^^ = ^^A(r)/87r (4.6) 

where p, and are the proper energy density, radial and tangential pressures 
respectively and r]i is the unit space-like vector on which the condition to be imposed 
is rjirj' = — 1. Here Pr is the pressure in the direction of r]i whereas p±_ is the pressure 
on the two-space orthogonal to r]i. 
Now, for the spherically symmetric metric 

ds'' = e^'Ue - e^^'Ur'' - r\d0^ + sinHd(t>''), (4.7) 
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the Einstein-Maxwell field equations (4.1) - (4.6) corresponding to anisotropic charged 
fluid with spatially varying cosmological constant i.e. A = A(r), are given by 

e-^(A7r - l/r^) + l/r^ = 87rT° = Svrp + (4.8) 

e-\v'lT + 1/r^) - = -SttT^ = 87rp^ - , (4.9) 
e~\v"l2 + i/'V4 - i/'A74 + (z/' - A')/2r] 

= -SttT^ = -87rr^3 = SttpI + (4.10) 

and 

P^]' = 47^rW/^ (4.11) 

where the intensity of electric field, is defined as £^ = —e~^^^'^^l'^<^' and can 
equivalently be expressed, from equation (4.11), as 

E = \ r iTTr^ae^^^dr. (4.12) 

Here prime denotes differentiations with respect to the radial coordinate r only. 
In the above equations (4.8) - (4.10) we have considered that 

p = p + A(r)/87r, (4.13) 

Pr = Pr - A{r)/87r (4.14) 

and 

P± = P± - A(r)/87r, (4.15) 

where p, Pr and p± are the effective energy density, radial and tangential pressures 
respectively. 

The equation of continuity T*^; i = 0, is given by 
where q is the charge on the spherical system. 

We assume the relation between the radial and tangential pressures (Herrera and 
Varela 1994) as 

P±-Pr = Q;?V, (4.17) 
50 



where a is a constant. 

Hence, by use of equations (4.14) and (4.17), the equation (4.16) reduces to 

Now, equation (4.8) can be expressed in the following form as 

e-^ = 1 - 2M/r, (4.19) 
where the active gravitational mass, M, is given by 



M = 47r 







r^dr. (4.20) 



4.3 The solutions 
4.3.1 Model for p + = 

Let us now solve the equation (4.18) under the assumption between the stress-energy 
tensors as T^i = T°o, which implies that 

p + p^ = p + p^ = 0. (4.21) 

In order to make the equation (4.12) integrable we assume that 

(7 = (Joe-^/^ (4.22) 

where ctq is the charge density at r = of the spherical distribution, i.e. the central 
density of charge. 

Now, using condition (4.22) in equation (4.12), we get for the expression of 
electric charge and intensity of the electric field as 

4 

q = Er^ = -Traor^. (4.23) 

With the help of equations (4.21) and (4.23), the equation (4.18) reduces to 

^ = ^TTCToV + 2Q;gV, (4.24) 
ar 3 
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Thus the solution set is given by 

e-^ = =1- ^7rVoV2(5a2 - 2r'^) - ^naq^r'^{ha^ - Sr^), (4.25) 
45 15 

Pr = -{aq' + \i^al){a' - r') + (4.26) 

= aq\' - {aq' + ^7rao^)(a^ " + ^ (4-27) 



and 



.2 , 2 2\/„2 A(r) 



p(r) = (aq' + -na',){a' - r') - (4.28) 



The active gravitational mass which is defined in the equation (4.20), then, by virtue 
of the equations (4.23) and (4.28), takes the form as 

M(r) = -^7rVoV^[87rQ;a6(5a2 - Sr') + 3(50^ - 2r')]. (4.29) 
135 

Thus, the metric potentials A and u are given by 

= e- = 1 - H^. (4.30) 
r 

The total effective gravitational mass can be obtained, after smoothly matching of 
the interior solution to the exterior Reissner-Nordstrom solution on the boundary, 
as 

m = M{a) + ^ ^7rVo2a^(l + ^Trcva^), (4.31) 

2a 45 3 

which corresponds to the second case (B) of Herrera-Varela model (1994) and repre- 
sents "electromagnetic mass" model such that gravitational mass of a charged fiuid 
sphere originates from the electromagnetic field alone (Lorentz 1904; Feynman et 
al. 1964). This again corresponds to the Tiwari-Rao-Kanakamedala model (1984) 
with « = case and thus the present model reduces to isotropic one. 
Now, considering the observed vahics of mass, charge and radius of the electron (in 
relativistic units) as m = 6.76 x 10~^^ cm, q = 1.38 x 10~^^ cm and a = 10~^^ cm 
the value of a, from the equation (4.31), is given by 

a = -4.77 X 10^^ cm-^ (4.32) 
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For the above value of a, the energy density in equation (4.28) becomes 



p(r) = -6.81 X 10^^(0^ - r^) 



A(r) 



(4.33) 



The central energy density, po, at r = 0, then can be calculated as 



po = -6.81 X 10-^ - 



(4.34) 



Thus, from the equation (4.34) one can see that for Aq > the energy density of 
the electron is a negative quantity. It is to be noted here that in the cosmological 
context A positive is related to the repulsive pressure and hence an acceleration 
dominated universe as suggested by the SCP and HZT project report (Perlmutter 
et al 1998; Riess et al. 1998; Filippenko 2001; Kastor and Traschen 2002). However, 
equation (4.34) indicates that this negativity of energy density is also obtainable 
for Aq < (which indicates a collapsing situation of the universe (Cardenas et al. 
2002)) for its very small value. In this context it is also possible to show that at an 
early epoch of the universe when the numerical value of negative A is higher than 
that of the first term of p (i.e. ~ 10^^ at r = 0) obviously energy density is a positive 
quantity. Thus, in the case of decreasing negative value of A it is clear that there is 
a smooth crossover from positive energy density to a negative energy density via a 
phase of null energy density! However, these results confirm the vacuum equation 
of state p-\-pr — ^ (Davies 1984; Blome and Priester 1984; Hogan 1984; Kaiser and 
Stebbins 1984). 

We can also see that on the boundary, r — a, the total energy density becomes 



which shows its clear dependency on the cosmological constant. However, for A^ > 0, 
pa is negative whereas for A^ < 0, pa is as usual a positive quantity. Through a sim- 
ple and interesting exercise (as all the parameters related to the electron are known) 
one can find out the numerical value of A^, at the boundary of the spherical system 
from the equation (4.35), which equals ~ 10~^cm~^. This constant value of A^ is 
too large and might be related to an early epoch of the universe. Here for finding 



__K 

87r' 



(4.35) 
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out the total energy density pa it is considered that pa < Paverage, where Paverage is 
equal to m/|7ra^ as the energy density of the spherical distribution is decreasing 
from centre to boundary. 



4.3.2 Model for p + ^ 

Now using equation (4.23) the equation (4.18) can be written as 



d 
dr 



Pr - 



El 
Stt 



+ 2aq r. 



Assuming that the radial stress-energy tensor T\ = 0, one gets 



1) 



and 



Pr 



Stt' 



Using equations (4.37) and (4.38) in equation (4.36), we have 



^-1 



P + Pt= P + Pr^ 

Thus, equation (4.20) takes the form as 

^(4agV + EVtt 



M = 47r 



r dr. 



To make equation (4.40) integrable we assume that 

= 7rk{e^ - 1)(1 - R^) - 47ragV2, 

where A; is a constant and R = r/a, a being the radius of the sphere. 
Thus, the solution set is given by 

e-^ ^1- AR^(5-3R^), 

^{l- 2Af/^e^/''exp[5Btan-^B{6R^ - 5) - Itan'^B], 



(4.36) 



(4.37) 



(4.38) 



(4.39) 



(4.40) 



(4.41) 



(4.42) 
(4.43) 



54 



and 



Pr = \Ke^ - 1)(1 - R') - \aqh' + (4.44) 
= \k{e>^ - 1)(1 - R') + ^agV^ + (4.45) 



p = fc(l - i?2)[l - \{e' - 1)] + iagV - (4.46) 



where the constant A = 8nka^/15. 

By apphcation of the matching condition at the boundary we again get the total 
effective gravitational mass, which in the present case takes the form 

8 0^ 

m = —nka^ + — . (4.47) 
15 2a ^ ^ 

In view of the equation (4.41), for vanishing charge the constant k vanishes and 
hence makes the gravitational mass in the equation (4.47) to vanish. Thus, the 

present case p + pr = k{l — R^) ^ also represents electromagnetic mass model. 
Now, the constant k can be expressed in terms of the known values of the electric 
mass, radius and charge as 

At r = 0, the energy density, from the equation (4.46), is then given by 

po = -5.68 X 10"^ - ^. (4.49) 

As, in the case of electron, /c is a negative quantity so for Aq > the central 
energy density po is negative only. However, for Aq < the central energy density 
may respectively be negative and positive depending on the numerical value of k 
whether it is higher and lower than that of Aq. 
At r = a, the total energy density is given by 

Pa = -4.54 X 10-5 - ^. (4.50) 

OTT 

Similarly, for A^ > 0, the energy density is negative whereas for Aa < 0, it may 
either be negative or positive depending on the numerical value of Aa as discussed 
in the previous case. 
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4.3.3 A test model 



In the previous two cases we have quahtatively discussed the effect of cosmological 
parameter A(r) on the energy density p(r) of the electron. Let us now explore some 
quantitative effect and hence treat the equation (4.24) in a different way. If we 
substitute the value of p^; from equation (4.14), then integrating equation (4.24) we 
get 

Ae// = A(a) - A(r) = 87rp(r) - M{aq^ + ^Tra^) (a^ - r^) , (4.51) 

where Ag// is the effective cosmological parameter. 
We study the following cases: 

For the central value of the energy density of the spherical distribution, i.e. at r = 0, 
the effective cosmological parameter becomes 

Klff = A(a) - A(0) = Sttpo - 877(0;?^ + ^7r(7g)a^. (4.52) 

Considering that po ^ P average the effective cosmological parameter, at r = 0, for 
the proper numerical values of the charge and radius of the electron can be found 
out as 

Klff = 1.71 X lO-^cm-2. (4.53) 

On the other hand, at the boundary, r = a, of the spherical distribution the effective 
cosmological parameter becomes 

A:^^ = A(a)-A(a) = 0. (4.54) 

Thus, from the equations (4.53) and (4.54) it is shown that the effective cosmological 
parameter has a finite value at the centre of the electron which decreases radially 
and becomes zero at the boundary. 

4.4 Discussions 

We see from the above analysis that the cosmological parameter A has a definite 
role even on the energy density of micro-particle, like electron. We, therefore, feel 
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that it may also be possible to extrapolate the present investigation to the massive 
astrophysical bodies to see the effect of spatially varying cosmological parameter on 
their energy densities and vice versa. 

The proper pressure Pr, in general, being positive as evident from the equation 
(4.26) is in accordance with the condition (4.21) which may be explained as due 
to vacuum polarization (Gr0n 1985). In this connection it is mentioned by Bonnor 
and Cooperstock (1989) that the negativity of the active gravitational mass and 
hence negative energy density for electron of radius a ~ 10~^^ is consistent with the 
Reissner-Nordstrom repulsion. We would also like to mention here that the equation 
of state in the form p+p = is discussed by Gliner (1966) in his study of the algebraic 
properties of the energy-momentum tensor of ordinary matter through the metric 
tensors and called it the p-vacuum state of matter. It is also to be noted that the 
gravitational effect of the zero-point energies of particles and electromagnetic fields 
are real and measurable, as in the Casimir Effect (1948). According to Peebles and 
Ratra (2002), hke all energy, this zero-point energy has to contribute to the source 
term in Einstein's gravitational field equation. This, therefore, demands inclusion 
of vacuum energy related term cosmological constant in the field equation. In this 
regard it is interesting to recall the comment made by Einstein (1919) where he 
stated that "... of the energy constituting matter three-quarters is to be ascribed 
to the electromagnetic field, and one-quarter to the gravitational field" and did 
"disregard" the cosmological constant in his field equation is in contradiction to the 
present result as shown in the equation (4.31) and (4.47). 



The contents of this chapter published in Physics Letters A (2004) 322 150. 
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Chapter 5 



Relativistic Electromagnetic Mass 
Models: Charged Dust 
Distribution in Higher Dimensions 



"The most incomprehensible thing about 
the world is that it is comprehensible. " 
- Albert Einstein 



5.1 Introduction 

In many theories higher dimensions play an important role, specially in super string 
theory (Schwarz 1985; Weinberg 1986) which demands more than usual four dimen- 
sional space-time. This is also true in studying the models regarding unification of 
gravitational force with other fundamental forces in nature. In the case of a simple 
solution to the vacuum field equations of general relativity in (4 -|- 1) space-time 
dimensions Chodos and Detweiler (1980) have shown that it leads to a cosmology 
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which at the present epoch has (3 + 1) observable dimensions in which the Einstein- 
Maxwell equations are obeyed. Lorenz-Petzold (1984) has studied a class of higher 
dimensional Bianchi-Kantowski-Sachs space-times of the Kaluza-Klein type whereas 
Ibanez and Verdaguer (1986) have examined radiative isotropic cosmologies with ex- 
tra dimensions related to FRW models. 

In this connection it is interesting to note that electromagnetic mass models where all 
the physical parameters, including the gravitational mass, are arising from the elec- 
tromagnetic field alone have been extensively studied (Tiwari et al. 1984; Gautreau 
1985; Gr0n 1986; Ponce de Leon 1987; Tiwari and Ray 1991a; Ray and Bhadra 
2004a,b) in the space-time of four dimensional general relativity. Thus it is believed 
that study of electromagnetic mass models in higher dimensional theory will be 
physically more interesting. 

Under this motivation we have considered here a static spherically symmetric charged 
dust distribution corresponding to higher dimensional theory of general relativity. 
It is proved, as a particular case, from the coupled Einstein-Maxwell field equations 
that a bounded and regular interior static spherically symmetric charged dust so- 
lution, if exists, can only be of purely electromagnetic origin. An example, which 
is already available, is examined in this context and is shown that the solution set 
satisfies the condition of being electromagnetic origin. 

5.2 The Einstein-Maxwell Field Equations 

The Einstein-Maxwell field equations for the case of charged dust distribution are 
given by 



(5.1) 



1/2 



(5.2) 



and 




(5.3) 
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where F^^ is the electromagnetic field tensor and J* the current four vector which is 
equivalent to 



(5.4) 



(7 being the charge density and is the four velocity of the matter satisfying the 
relation 



UjU^ — 1. 



The matter and electromagnetic energy momentum tensors are given by 



rpl 



pU Uj 



and 



rpi (em) 1 r „ „i TP pWl 



where p is the proper energy density. 

Now we consider the {n + 2) dimensional spherically symmetric metric 

ds2 = e'^W(it2-e^M(ir2-r2X„2 



where 

The convention adopted here for coordinates are 



n-l 



x^^r, x^^9i, x^^92, = = i 



and also 



A 2 2-2/1 

9u = -e , 5^22 = -r , 5^33 = -r sin 9i, 

"n— 1 

,2 •_2/) -2/ 



5^44 = sin 6'ism 02, g(n+i){n+i) = 



(5.5) 



(5.6) 



(5.7) 



(5.8) 



(5.9) 



5'(n+2)(n+2) = 6 



(5.10) 
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As we have considered here a static fluid, so 

u' = [0,0,0, (n + l)times, e'"/^] 



(5.11) 



and 

^J^^J^^ ...r+^ = 0, J"+2 ^ (5.12) 

so that the only non- vanishing components of of equations 5.2 and 5.3 are i^i(n+2) 
and F(„+2)i- 

In view of the above, the Einstein-Maxwell field equations for static spherically 
symmetric charged dust corresponding to the metric (5.8) are 

e-\nu'/2r + n{n - l)/2r'^] - n{n - l)/2r^ = -E^, (5.13) 

e-\u"/2 + u'^A - iy'X'/4 + (n - l)(z/' - A')/2r + 

{n - l){n - 2)/2r'] - {n - l){n - 2)/2r^ = E^, (5.14) 

e~^[nX'/2r-n{n-l)/2r^]+n{n-l)/2r^ = 8np + E^ (5.15) 

and 

[r^'E]' = 47rr"ae-^/2 ^gjg^ 

where E, the electric field strength, is defined as = —e'^^^'^^^/'^ct)' , the electrostatic 
potential being related to the electromagnetic field tensor as F(n+2)i = — Fi(„+2) = 

CP'. 

5.3 Higher Dimensional Electromagnetic Mass Mod- 
els 

^From the field equations (5.13) and (5.15), we have 

e-^{v' + A') = 167rrp/n. (5.17) 
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Now, from the conservation equations T^j;i — 

pu' = [q^^/Anr^ + (n - 



one obtains 

2)E727rr 



(5.18) 



where the charge, g, is related with the electric field strength, through the integral 
form of the Maxwell's equation (5.16), which can be written as 

q^Er"^ 47r T (7e^/V"dr. (5.19) 
Jo 

Again, equation (5.15) can be expressed in the following form as 

e"^ = 1 - 4M/nr"-^ (5.20) 

where the active gravitational mass, M, is given by 

M^An r[p + E^/8n]r''dr. (5.21) 
Jo 

Hence, following the technique of Tiwari and Ray (1991a) we see that for the van- 
ishing charge density, a, of the equation (5.19) one can arrive, through the equation 
(5.18), at the unique relation 

pu' = 0. (5.22) 

Thus, we have the following two cases: 
Case I: p^O, i/' = 

For this case, from equations (5.13) and (5.14), we have A as a constant quantity. 
This in turn makes p equal to zero and hence by virtue of equation (5.22) space-time 
becomes flat. 

Case II: p = 0, u' ^ 
In this case, from equations (5.20) and (5.21), A becomes zero. Then from equation 
(5.17), we have the metric potential as a constant and again the space-time becomes 
fiat. 

We are not considering here the third case, viz. p — v' — 0, which is quite a trivial 
one. However, from the above cases (I) and (H) it is evident that, at least at a 
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particular case, all the charged dust models are of electromagnetic origin, viz., all 
the physical parameters are originating purely from electromagnetic field. This type 
of models are known as electromagnetic mass models in the literature (Lorentz 1904; 
Feynman et al. 1964). 



An example: 



The solution set obtained by Khadekar et al. (2001) for the static spherically 
symmetric charged dust is as follows: 



(5.23) 



in -I] 



and 



a 



N+{n-l) 

_ Nn{n - If 
~ 87rr2[7V+ (n- 1)]2 



(5.24) 



(5.25) 



(5.26) 



47r2VV2[Ar + (n- 1)]2" 

where A and N both are constants with the restriction that N >0. 

The total charge and mass of the sphere in terms of its radius, a, are respectively 

given by 

N[n{n- 1)^/^0.''-'^ 



Q 



2V2[Ar+ (n- 1) 



and 



Na 



n-l 



m — 



N + {n-l)' 

The charge and mass densities in the present case take the relationship 



(5.27) 



(5.28) 



a = [2(1 - l/n)]VV 



(5.29) 
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Therefore, the charge and mass densities are proportional to each other with the 
constant of proportionality [2(1 — which takes the value unity for n = 2 i.e. 

in the four dimensional case and the relation (5.29) reduces to the usual form 

a = ±p (5.30) 

which is known as the De-Raychaudhuri (1968) condition for equilibrium of a charged 
dust fluid. 

Thus, from equations (5.23) - (5.29), it is evident that all the physical quantities 
including the effective gravitational mass vanish and also the spherically symmetric 
space-time becomes flat when the charge density vanishes implying N = 0. The so- 
lution here, therefore, satisfies the criteria of being of purely electromagnetic origin. 

5.4 Discussions 

The present paper is, in general, higher dimensional analogue of the work of Tiwari 
and Ray (1991a) whereas the example given here (Khadekar et al. 2001) is the 
higher dimensional analogue of the paper of Pant and Sah (1979). Thus, we have 
presented here a model which corresponds to spherically symmetric gravitational 
sources of purely electromagnetic origin in the space-time of higher dimensional 
theory of general relativity. It has been already proved in the four dimensional 
case of the present paper (Tiwari and Ray 1991a) that a bounded continuous static 
spherically symmetric charged dust solution, if exists, can only be of electromagnetic 
origin. Hence, this is also true in the higher dimensional case in general theory of 
relativity. 

In this regard we would like to discuss briefly the role of higher dimensions in 
different context. It have been shown by Ibanez and Verdaguer (1986) that for 
the open models related to FRW cosmologies the extra dimensions contract as a 
result of cosmological evolution whereas for flat and closed models they contract 
only when there is one extra dimension. Fukui (1987) recovers Chodos-Detweiler 
(1980) type solutions, as mentioned in the introduction, where the Universe expands 
as i^/^ by the percolation of radiation into 4D space-time from the flfth dimension. 
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mass, although the 5D space-time-mass Universe itself is in vacuum as a whole. It is 
interesting to note that considering mass as fifth dimension a lot of other works also 
have been done by several researchers (Wesson 1983; Banerjee, Bhui and Chatterjee 
1990; Chatterjee and Bhui 1990) which contain Einstein's theory embedded within 
it. In one of such investigations it is argued that a huge amount of entropy can 
be produced following shrinkage of extra-dimension which may account for the very 
large value of entropy per baryon observed in 4D world (Chatterjee and Bhui 1990). 
Kaluza-Klein type higher dimensional inflationary scenario have been discussed by 
Ishihara (1984) and Gegenberg and Das (1985) where it is shown that the contraction 
of the internal space causes the inflation of the usual space. 



The contents of this chapter published in Astrophysics and Space Science (2006) 
302 153. 
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Chapter 6 



Relativistic Anisotropic Charged 
Fluid Spheres with Varying 
Cosmological Constant 



"What are man's truths ultimately? 
Merely his irrefutable errors!" 
- Nietzsche 



6.1 Introduction 

The cosmological constant A, related to the energy of space, introduced by Einstein 
in general relativity has become very significant from the view point of cosmology. 
Though Einstein ultimately abandoned it stating that it was a "blunder" in his life 
but Tolman keeps it as a constant quantity in his field equations even in 1939's fa- 
mous work related to astrophysical system. It is also, in favor of keeping A, argued 
by Peebles and Ratra (2002) that like all energy, the zero-point energy related to 
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space has to contribute to the source term in Einstein's gravitational field equa- 
tions. However, it is being gradually felt that the erstwhile cosmological constant A 
is indeed a scalar variable dependent on time rather than a constant as was being 
believed earlier. Recently, this variation in cosmological constant is also observation- 
ally confirmed due to the evidence of high redshift Type la supernova (Perlmutter 
et al. 1998; Riess et al. 1998) for a small decreasing value which is < 10~^^ cm~^ 
at the present epoch. Obviously, once the A becomes a scalar its dependence need 
not be limited only to time coordinate (as in cosmology). Since it enters in the 
field equations as a variable, it must be dependent on space coordinates as well. 
Therefore, in general, the A is a scalar variable depending on, either or both, space 
and time coordinates. It is argued that just as in cosmology the dependence of A 
on time has been found to be of vital importance playing a significant role now, 
its dependence on space coordinates is equally important for astrophysical problems 
(Chen and Wu 1990; Narlikar, Pecker and Vigier 1991; Ray and Ray 1993; Tiwari 
and Ray 1996). 

With this view point, we consider here an anisotropic charged static fiuid sphere 
by introducing a scalar variable A dependent only on the radial coordinate r. The 
field equations thus obtained, under certain mathematical assumptions, yield a set 
of solutions which has another historical importance, known in the name of Electro- 
magnetic Mass Models (EMMM) in the literature (Lorentz 1904; Hoffmann 1935; 
Feynman, Leighton and Sands 1964; Tiwari, Rao and Kanakamedala 1986; Wilczek 
1999). The cfi^ectivc gravitational mass of these models depends on the electro- 
magnetic field alone, viz., the effective gravitational mass vanishes when the charge 
density vanishes. Such models have been studied by several authors (Ray and Ray 
1993; Tiwari, Rao and Kanakamedala 1984; Tiwari, Rao and Ray 1991; Gautreau 
1985; Gr0n 1985, 1986a,b; Ponce de Leon 1987a,b, 1988; Tiwari and Ray 1991a,b, 
1997; Ray, Ray and Tiwari 1993). All these EMMMs, however, have been obtained 
under a special assumption p + p — 0, where p is the matter-energy density and p 
is the fiuid pressure under the general condition that p > and p < 0. This type of 
equation of state, viz., p — ^p with 7 = — 1, implies that the matter distribution is 
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in tension and hence the matter is known, in the hterature, as a 'false vacuum' or 
'degenerate vacuum' or 'p-vacuum' (Davies 1984; Blome and Priester 1984; Hogan 
1984; Kaiser and Stebbins 1984). A natural question arises whether there exists any 
EMMM where this condition is violated, i.e., when p + p 7^ 0. This is the main mo- 
tivation of the present investigation and here we have shown that even ion: p + p ^ 
EMMM can be constructed. However, the same question was addressed by Tiwari, 
Rao and Ray (1991) and obtained EMMM in the isotropic and axially-symmetric 
matter distribution for charged dust case only whereas Ray and Bhadra (2004b) 
searched for solution to the problem by employing a relation between the radial and 
tangential pressures as p± = Pr + aq^r"^ where a is a non-zero constant factor and q 
is electric charge of the spherical system of radius r. We shall consider in the present 
investigation different relation and would like to observe that how this helps us to 
find out several class of solutions related to EMMM. 



6.2 The Einstein-Maxwell Field Equations 

Let us consider a spherically symmetric line element 

ds' = Qijdx'dx^ = e^^'Mt" - e^^^'^dr^ - r\d9^ + sin^9d(j>^), = 0, 1, 2, 3) (6.1) 

where u and A are the metric potentials. 

Now, the Einstein field equations for the case of charged anisotropic source are 

G'j = R'j - ^g'jR = -K[r/""^ + + T V^], (6.2) 

where T^^J^\ T^^^^^ and T"*^"""^^ are respectively the energy- momentum tensor com- 
ponents for the anisotropic matter source, electromagnetic field and vacuum. The 
explicit forms of these tensors are given by 

T't^ = (p + - p^g'i + (p± - Pr)v%, (6.3) 
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and 

~ Stt 

with UiU^ — —riirf — 1. The Maxwell electromagnetic field equations are given by 



T'T"' = (6-5) 



[{-9f"F%^ = A^J\-gf'' (6.6) 

and 

FiijM = 0' (6-7) 
where the electromagnetic field tensor Fij is related to the electromagnetic potentials 
through Fij = Aij — Aj^i which, obviously, is equivalent to the equation (6.7). 
Further, is the 4- velocity of a fiuid element, J' is the 4-current satisfying J' = au^, 
where a is the charge density, and k — Sir {in relativistic unit G = C = 1). Here 
and in what follows a comma denotes the partial derivative with respect to the 
coordinates (involving the index). 

The Einstein-Maxwell field equations (6.2) - (6.7) corresponding to static anisotropic 
charged source with cosmological variable, are then given by 

e-^{X'/r - l/r^) + l/r^ = 87rr% = Snp + + A, (6.8) 
e'^{u'/r + l/r^) - l/r^ = -SttT^ = Snpr - - A, (6.9) 
e-^[u"/2+u'^/4-u'X'/4+{u'-X')/2r] = -SttT^s = -SttT^ = Snp^+E^-A (6.10) 
and 

(r^£;)' = 47rrW/^ (6.11) 

The equation (6.11) can equivalently, in terms of the electric charge q, be expressed 
as 

q(r) = r^E(r) = /' inr^ae^/^dr (6.12) 
Jo 

where Pr, p± and E are the matter-energy density, radial and tangential pressures 
and intensity of the electric field respectively. Here prime denotes derivative with 
respect to radial coordinate r only. 
The equation of continuity, T^j; i = is given by 

4- \pr - {E^ + A(r)}/87rl + {p + Pry/2 = Ey27rr + 2{p^ - p,)/r. (6.13) 
ar L J 
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Now, our spherically symmetric fluid distribution under investigation is anisotropic 
in nature. This means that the radial and tangential pressures are, in general, 
unequal so that the simplest relation between them we assumed here as 

p± = npr, (ny^l). (6.14) 

Assuming further that the radial stress T^i — (Florides 1977; Kofinti 1985; Gr0n 
1986a; Ponce de Leon 1987b; Tiwari and Ray 1996) one gets 

i/'=(e^-l)/r (6.15) 

and 

Pr = [E^ + A{r)]/87r. (6.16) 
Using equations (6.14) - (6.16), in equation (6.13), we get 

p + p, = [(n+l)E' + (n-l)A(r)]/27r(e^-l). (6.17) 

Similarly, equations (6.8) and (6.9) yield 

e-^{u' + A') = 87rr(p + Pr). (6.18) 

Again, equation (6.8) through the equation (6.16), gives 

e-^ = 1 -2m(r)/r, (6.19) 

where m(r), called the efi^ective gravitational mass, takes the form 

rr 

mir) = Mir) + jiir) = A-k / [p + Pr]r^dr, (6.20) 

Jo 

the active gravitational mass of Schwarzschild type and the mass equivalence of 
electromagnetic field, respectively, being defined as 

M(r) = 47r / pr'^dr, i^{r) = Att f p^dr. (6.21) 
Jo Jo 

Now, from the above equation (6.20) it is easily observed that the condition p+Pr — 
0, yields a flat space-time through the equations (6.19) and (6.15) and has been 
considered by Tiwari, Ray and Bhadra (2000) in another context. Hence, the non- 
trivial solutions exits here for the case p + Pr ^ only. 
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6.3 Solutions for the Static Charged Fluid Spheres 



Now we solve the equation (6.17) under the constraint p + 7^ 0, assuming differ- 
ent mathematical conditions. As the equation (6.17) is involved with two physical 
parameters, E and A, so unless we specify one parameter in terms of other it is not 
possible to solve the equation (6.20). Therefore, a plausible straightforward relation 
between these parameters may be of the form A(r) = ±£'(r)^. Physically this means 
that vacuum energy has contribution from the electrostatic field energy and proposal 
of this kind, in an implicit way, is not at all unavailable in the literature (Ray and 
Bhadra 2004b). Let us, therefore, assume the following two cases A(r) — E"^ — NAq 
and A(r) = —E^ + NAq which will yield solutions with physically interesting fea- 
tures as the analysis of the next section demonstrates it clearly. However, in both 
the cases we have taken the assumptions in a way so that the cosmological variable 
A does not vanish rather may be at most equal to Aq, the erstwhile cosmological 
constant having a finite non-zero value. Without considering Aq we will have A = 
at the boundary r = a which is a bit unphysical and may create difficulties, such as 
to entropy like problems (Beesham 1993). 

Case I 

Let us assume that the cosmological constant is a function of radial distance such 
that 

A(r) ^E^- NAo (6.22) 

where is a free parameter and Aq is the erstwhile non-zero cosmological constant. 
With the help of equations (6.17) and (6.22), the equation (6.20) takes the form 

m(r) = 2 r[2nE^ - (n - l)A^Ao]r^dr/(e^ - 1). (6.23) 
Jo 

To make equation (6.23) integrable we further assume that 

= q^/r^ = [fc2(e^ - 1)(1 - R^) + (n - l)A^Ao]/2n, (6.24) 

where A; is a constant and R — r/a, a being the radius of the sphere. This partic- 
ular choice for the electric intensity generates a model for charged sphere which is 
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physically very interesting as it is related to EMMM as will be seen later on. Thus, 
the solution set is given by 

e-^M = l-yli?2(5_3i?2)^ (6.25) 

e^W = (1 - 2A)^/^e^('')/^exp[5Ston-^S(6i?2 - 5) - ton-^S/2], (6.26) 

Pr{r) = p±{r)/n = [k'^{e^{r) - 1}(1 - R^) - NAo]/8nn (6.27) 

and 

p(r) = [k'^{l + 4n - e^(r)}(l - R^) + NAol/Snn, (6.28) 

where 

A = Ak^a^ld, B = [A/{12 - 2bA)Y'^. (6.29) 

Now, the exterior field of a spherically symmetric static charged fiuid distribution 
described by the metric (6.1) is the unique Reissner- Nordstrom solution given by 

- P{de'^ + sin'^ed(f)^). (6.30) 

Then, application of the matching condition on the boundary r — a yields the total 
effective gravitational mass in the following form 

m(a) = ^k'a' + (6.31) 
io Za 

Hence, in terms of the total gravitational mass, the total electric charge and the 
radius of the sphere, all the constants k, A and B can be expressed as 

k'^ = 15(2am - q^) /8a^ (6.32) 

A ^ {2am - q^) /2a^ (6.33) 

and 

= (2am - g^)/[24a2 - 25(2am - q% (6.34) 

Comparing the equation (6.20) (and subsequently via the equation (6.21)) with the 
above equation (6.31) we can easily recognize the first term (4/c^a^/15) in the right 
hand side as the total Schwarzschild mass whereas the second term [q^ /2a) is the 



ds''^ 



2m q 
1 + 



df 



2m q^ 
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total mass equivalence of the electromagnetic field. Now, k, in the first term of 
equation (6.31), is implicitly related to the charge q which is evident from the equa- 
tion (6.24). Therefore, in principle, the gravitational mass m expressed in equation 
(6.31) is of purely electromagnetic in origin. 

Case II 
Here our choice is 

Mr) = -E'^ + NKo. (6.35) 

Then, from the equations (6.17) and (6.35), the gravitational mass (equation (6.20)) 
reduces to 

mir) = 2 r\2E'' + (n - l)A^Ao]r'dr/(e^('^) - 1). (6.36) 
For the further assumption 

E"" = q^/r^ = [k^e^^^-^ - 1)(1 - R^) - {n - l)A^Ao]/2, (6.37) 
we obtained the solution set as 

e-^M = l-AR^{b- 3R^), (6.38) 

e^W = (1 - 2Af/^e^^'^/^exp[5Btan-'B{6R' - 5) - tan-'B/2\, (6.39) 

p{r) = NAo/8n (6.40) 

and 

p{r) = [4k\l - R^) - iVAo]/87r. (6.41) 

We see that A and retain the same form as in the case I and hence the total 
gravitational mass is also given by the equation (6.31). Further, it is observed that 
the present case automatically reduces to an isotropic one as the pressure p does 
not associated with the anisotropic factor n. 
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6.4 Physical Properties of the Static Charged Fluid 
Spheres 

Case I: K{r) ^ - NAq 

Subcase (1): The equation (6.24) indicates that for getting a direct dependence 
of k upon q one can admit the following relaxation such that {\) N = Q when Aq 7^ 

and (n — 1) 7^ 0, (ii) (n — 1) = when 7^ and Aq 7^ and (iii) Aq = when 
(n — 1) 7^ and A^ 7^ 0. The third possibility seems to contradict the observational 
results related to the Supernova type la where the cosmological constant is found 
to be a non-zero positive value whereas the second one provides an isotropic case. 
Then, suitably opting for A^ = one obtains 

q{r f = k^r\e^('^ - - R^) /2n, (6.42) 

Pr{r) = P±{r)/n = A;'(e^(") - 1)(1 - R')/Snn (6.43) 

and 

p(r) - k'^{l + 4n - e^('')}(l - R^)/87m. (6.44) 

Thus, for vanishing electric charge the gravitational mass in the equation (6.31), 
including all the physical parameters (viz., pressures and density), vanishes and one 
obtains EMMM. 

Subcase (2): The central and the boundary pressures are found to be equal here 
i.e., Pr{0) = p±{0)/n = —NAo/Snn and Pr{a) = p±{a)/n = —NAo/87rn respec- 
tively whereas the respective densities are (4n/c^ -|- NAq) /Sim and NAo/Snn. So, 
the present model has a constant pressure throughout the sphere though the density 
decreases from centre to boundary. For the value N — 0, however, we have zero 
pressure, both at the centre and boundary. The density decreases from the non-zero 
central value A;^/27r and then smoothly decreases to zero at the boundary. Thus, 
with N — the present model goes to a physically well-behaved static charged dust 
case. 
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Subcase (3): In the above analysis for the general value of we have seen that the 
fluid pressure is altogether negative whereas the density is a positive quantity. Now, 
it can be observed from the equations (6.27) and (6.28) that p -\- Pr ^ Q except at 
the boundary where it is equal to Pr{o) — —p{o) such that p{a) > and p{a) < 
as seen earher. Otherwise, it will have a general value A;^(l — i?^)/27r . The central 
value is then ^,.(0) = — p(0) where p(0) = p(0) — A;^/27r. As the model demands for 
p > and p < 0, so the condition to be satisfied here is p(0) > A;^/27r. The general 
condition for the negative pressure and positive density is then p > k^{l — i?^)/27r 
for all r < a. These results are also true for the sub-case N — 0. 



Case II: A(r) = -E"^ + NAo 

Subcase (1): Here iov N — Q the solution set, as obtained in the equations (6.38), 
(6.40) and (6.41), reduces to 

p(r) = (6.45) 

and 

p(r) = A;2(l-i?2)/27r, (6.46) 
when the electric charge is given by 

q{rf = A;2(e^('-) - 1)(1 - R^y /2. (6.47) 

Thus, as in the previous case, for g = wc get k = which in turn makes mass, 
pressure and density to vanish and also the space-time becomes flat. Thus, the 
model presented here is an EMMM. 

Subcase (2): In the present case also the central and the boundary pressures are 
equal with a value NAq/Stt and the respective densities are {Ak"^ — NAo)/8'k and 
—NAq/Stt. The pressures become zero for N — 0, a,t the centre and boundary, and 
densities have the positive central value /c^/27r whereas the boundary value is zero. 
Thus, we again get a physically interesting charged dust case with N — which now 
corresponds to a case of isotropic fluid sphere. Here also the behavior is regular and 
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well defined. 

Subcase (3): In the present case also, by virtue of equations (6.40) and (6.41), p+Pr 7^ 
which reads here os p -\- Pr — /c^(l — -R^)/27r. Here the central value, at r = 0, is 
p — —p where p — (p — A;^/27r) and the boundary one is p = —p. Due to negative 
value of the density here the condition on the pressure to be imposed is p > /c^/27r. 
Thus, the present situation, viz. Case II, clearly provides an EMMM even with a 
positive pressure and therefore contradicts the comment made by Ivanov (2002) that 
"... electromagnetic mass models all seem to have negative pressure." The same 
result, i.e. the positivity of pressures are also available in some cases of the work 
done by Ray and Das (2002) related to EMMM. However, the explanation given 
here is valid for any positive value of and so the situation could completely be 
opposite if one assigns on N any negative value. At this stage, we should not put 
any restriction on the choice of the value of A^. This is because, in general, for a 
fluid sphere we should have p > and p > so that the weak energy conditions are 
satisfied. But there are also some special situations available within the spherical 
system (particularly in the case of electron with the radius ~ 10~^^) where the 
energy condition is violated due to negative energy density (Cooperstock and Rosen 
1989; Bonnor and Cooperstock 1989; Ray and Bhadra 2004b). Thus, choosing the 
proper signatures of A", we can have a class of models with diverse characters. 



6.5 Role of A: Previous and Present Status 

The cosmological constant was introduced by Einstein in his field equation to obtain 
a static cosmological solution because of the fact that due to gravitational pull 
everything will collapse to a point and hence a un-wanting situation of singularity 
will take place. However, he was not satisfied with this new physical quantity as it 
seemed to violate Machian principle which he tried to incorporate in the framework 
of his General Relativity. He thus, ultimately rejected it mainly for two reasons: 
(i) that the theoretical work of de Sitter showing that the Einstein's field equations 
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admitted a solution for empty Universe and (ii) that the experimental discovery of 
expanding Universe by Hubble. 

As stated in the introduction, the concept of cosmological constant has been revived 
recently in the case of early Universe scenario and even in particle physics. It is 
gradually being felt that A, the erstwhile cosmological constant is available rather 
than a constant, as was being believed earlier, varying with space or time or both 
(Tiwari and Ray 1996; Ray, Ray and Tiwari 1993; Tiwari, Ray and Bhadra 2000). 
Further, this varying A may be positive or negative (by imposing the condition that 
its value is not equal to zero). For instance, according to Zel'dovich (1968) the 
effective gravitational mass density of the polarized vacuum is negative. Similarly, 
the equation of state p + p = 0, employed by Tiwari, Rao and Kanakamedala (1984) 
to construct EMMM as a solution of Einstein-Maxwell field equations, provides 
negative pressure. It may be emphasized here that positive density has significant, 
rather major role in inflationary cosmology whereas negative density has influence 
on elementary particle models. The gravitational mass inside the spherical charged 
body is negative for r < 5a/ A, where r is the radial coordinate and a is the radius 
of the sphere. It is argued by Gr0n (1986a,b) that this negative mass and the 
associated gravitational repulsion is due to the strain of the vacuum because of 
vacuum polarization. He also argued that if a vacuum has a vanishing energy, then 
its gravitational mass will be negative and the observed expansion of the universe 
may be explained as a result of repulsive gravitation. Now, if we consider a negative 
A having a repulsive nature as was considered by Einstein then this gets the same 
status of negative pressure and also can be identified with the Poincarc stress. This 
repulsive gravitation associated with negative A can also be explained as the source 
of gravitational blue shift (Gr0n, 1986a). On the contrary, positive A will be related 
to gravitational red shift. It may also be pointed out that according to Ipser and 
Sikivie (1984) domain walls are sources of repulsive gravitation and a spherical 
domain wall will collapse. To overcome this situation the charged "bubbles" with 
negative mass keep the wall static and hence in equilibrium. In this regard, we may 
also add that A, via repulsive gravitation, is related to domain walls and playing an 
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important physical role. 

Very recent observations conducted by the SCP and HZT (Perlmutter et al. 1998; 
Riess et al. 1998; Fihppenko 2001; Kastor and Traschen 2002) show that the present 
value of A is positive one and hence related to the repulsive pressure. It is believed 
that the present state of acceleration dominated Universe is due to the driven force 
of this A. It is, therefore, to be noted that the negative A corresponds to a collapsing 
situation of the Universe (Cardenas et al. 2002). 



6.6 Conclusions 

(i) In both the above cases / and //, it is possible to show that EMMM can be ob- 
tained, in principle, using the constraint p + p 7^ 0. This particular point remained 
unnoticed by Gr0n (1986a,b) and Ponce de Leon (1987a,b) both. 

(ii) It can be noted that in terms of energy-momentum tensor of the fluid the con- 
dition p + p = implies T^i — T^o^^ whereas p + p constraint may be expressed 
as T-*^! = as we have adopted in the present approach. It is also interesting to 
note that p + p = and hence T^i = T^o can be expressed in terms of the metric 
tensors (vide equation (6.1)) as 5'oofi'ii = ~1- A coordinate-independent statement 
of this relation is obtained by Tiwari, Rao and Kanakamedala 1984) by using the 
eigen values of the Einstein tensor C^j. 

(iii) We would like to mention here that the solutions obtained by Gr0n (1986a,b) 
and Ponce de Leon (1987a,b) represent a neutral system, viz., though the net charge 
is not zero but the charge on the surface of the spherical system vanishes. The mod- 
els of the present paper, in general, do not correspond to this situation because of 
the fact that the electric field and hence cosmological constant does not vanish at the 
boundary. In both the Case I and Case II, the values of electric field, respectively, 
are {n — l)NAo/2n and — (n — l)iVAo/2 whereas those for cosmological parameters 
are —{n+l)NAo/2n and (n + l)iVAo/2. Therefore, the present solutions correspond 
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to a charged fluid sphere. Of course, for N = 0, like Gr0n (1986a,b) and Ponce de 
Leon (1987a,b), we have neutral spheres (equation (6.42) of the Case / and equation 
(6.47) of the Case IT). Thus, we have a class of solutions related to charged as well 
as neutral systems depending on the values of N . 



The contents of this chapter has communicated to journal for publication. 
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Chapter 7 



Conclusions 



"So we come hack again to the original idea of Lorentz - may be all the 
mass of an electron is purely electromagnetic, may be the whole 0.511 
Mev is due to electrodynamics." 
- Feynman et al. (1964) 



Electromagnetic mass models which are the sources of purely electromagnetic 
origin "have not only heuristic flavor associated with the conjecture of Lorentz but 
even a physics having unconventional yet novel features characterizing their own 
contributions independent of the rest of the physics" (Tiwari 2001). This is, as Ti- 
wari (2001) guess "may be due to the subtle nature of the mass of the source (being 
dependent on the electromagnetic field alone)". Therefore, in our whole attempt we 
have tried to explore "the subtle nature of the mass of the source" . However, to 
do this under the general relativistic framework, we have considered Einstein field 
equations in its general form, i.e., with cosmological constant A which also acts as 
a source term to the energy-momentum tensor. If we consider that A has a variable 
structure which is dependent on the radial coordinate of the spherical distribution, 
viz., A = A(r) then it can be shown that A is related to pressure and matter energy 



80 



density. Hence it contributes to the effective gravitational mass of the system. 
It is seen that equation of state has an important role in connection to electromag- 
netic mass model. Therefore, at first we have obtained electromagnetic mass model 
under the condition p + p = 0. However, later on it is shown that electromagnetic 
mass model can also be obtained by using more general condition p-\-p ^Q. 
The model considered in our work, in general, corresponds to a charged sphere with 
cosmological parameter in such a way that it does not vanish at the boundary. The 
idea behind is that the cosmological parameter is related to the zero point vacuum 
energy it should have some finite non-zero value even at the surface of the bounding 
system. For this type of spherical system we can have a class of solutions related to 
charged as well as neutral configurations. 

It can be shown that these models have positive energy densities everywhere. Their 
corresponding radii are always much larger than 10~^^ cm. Furthermore, as the 
radii of these models shrink to zero, their total gravitational mass becomes infinite. 
It have been shown by Bonnor and Cooperstock (1989) that an electron must have 
a negative energy distribution (at least for some values of the radial coordinate). 
In this connection we have shown that the cosmological parameter A has a definite 
role on the energy density of micro particle, like electron. At an early epoch of 
the universe when the numerical value of negative A was higher than that of the 
energy density p, the later quantity became a positive one. In the case of decreasing 
negative value of A there was a smooth crossover from positive energy density to a 
negative energy density. 

So far we have referred electron to be a spherically symmetric distribution of matter 
deprived of spin and magnetic moment. As an alternative way both Bonnor and 
Cooperstock (1989) as well as Herrera and Varela (1994) suggest that both spin and 
magnetic moment can be introduced at classical level through the Kerr-Newman 
metric. However in this context it is to be mentioned here that the Kerr-Newman 
metric cannot be valid for distance scales of the radius of a subatomic particle. We, 
therefore, thought that the problem can be tackled in the frame work of Einstein- 
Cartan theory where torsion and spin are inherently present. In this case, the only 
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way is to take the spin to be the 'intrinsic angular momentum' that is the spin of 
quantum mechanical origin. In our work considering the spins of all the individual 
particles are assumed to be oriented along the radial axis of the spherical systems 
we have obtained some interesting solutions with physical validity. However, though 
our present approach via Einstein-Cartan theory to inject spin may be interesting 
it, at once, demands some alternative means to provide spin and magnetic moment. 
This may be possible through Dirac-Maxwell theory where spin and magnetic mo- 
ment are naturally incorporated through the Dirac spin. We would like to pursue 
this problem in future investigations. 

Another important point we would like to mention here that in all the previous inves- 
tigations wc have studied electromagnetic mass models in 4-dimensional Einstein- 
Maxwell space-times only. Therefore, one can ask whether electromagnetic mass 
models also can exist in higher dimensional theory of General Relativity. We have 
presented a model which corresponds to spherically symmetric gravitational sources 
of purely electromagnetic origin in the space-time of {n + 2) dimensional theory of 
general relativity. 

We have also taken up the problem of anisotropic fluid sphere as studied earlier in 
a different view point. By expressing A in terms of electric field strength E we have 
explored some possibilities to construct electromagnetic mass models using the con- 
straint p + p ^ 0. We would like to mention here that unhke the solutions of Gr0n 
(1986a,b) and Ponce de Leon (1987a,b) in the present investigation, in general, the 
electric field (and hence the cosmological constant) does not vanish at the bound- 
ary. However, it is shown that the class of solutions obtained here arc related to 
charged as well as neutral systems of Gr0n (1986a,b) and Ponce de Leon (1987a,b) 
depending on the values of the parameter A^. 

It is to be mentioned here that other than Dirac-Maxwell theory where spin and 
magnetic moment are naturally incorporated through the Dirac spin, some other pos- 
sibilities are awaiting to be investigated under the scheme of electromagnetic mass 
models. One of such possibilities is to study the relationship between the structures 
of sohton which have been identified with the electromagnetic field to that of elec- 
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trons which are also identified with the electromagnetic field (Tiwari 2001). This 
can be done by the use of Zakharov-Belinsky method to solve the Einstein- Maxwell 
equations. Another possibility is to conjecture that Weyl line-mass solutions and 
cosmic strings are identical entities, because it has been shown by Linet (1985) and 
Hiscock (1985) that the Weyl fine- mass solutions can be identified with the cosmic 
strings. On the other hand Weyl fine-mass solutions have been identified with the 
electromagnetic mass models (Tiwari et al. 1991). 
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